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ABSTRACT 


Using digitally processed MSS/LANDSAT data as auxiliary 
variable , a methodology to estimate wheat ( Triticum aestivum L) area 
by means of sampling techniques was developed. To per form this 
research , aerial photographs covering 720 km 2 in Crus Alta test site at 
the NW of Rio Grande do Sul State, , were visually analysed. LAUDS AT 
digital data obtained on Sept 2, 1979 were analysed using 
non-supervised and supervised classification algorithms j as 
post-processing the classification was submitted to the spactial 
filtering, To estimate wheat area, the regression estimation method was 
applied and different sample sizes and various sampling units (10, 20, 

30, 40 and SO km 2 ) were tested. Based on the four decision criteria 
estabhished for this research, it was concluded that: (a) as the size 
of sampling units decreased the percentage of sampled area required to 
obtain similiar estimation performance also decreased; (l) the lowest 
percentage of the area sampled for wheat estimation with relatively 
high precision and accuracy through regression estimation was 13,90 % 
using 10 km 2 as the sampling unit with this approximately 86 %, time 
and cost were saved during interpretation of airphotographs ; and (a) 
wheat area estimation by direct expansion (using only aerial photographs) 
was less precise and accurate when compared to those obtained by means 
of regression estimation . 
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CAPITULO 1 


introducao 


A necessldade de suprir a populagao de um paTs ondehares^ 
trigoes de area, clima e tecnologia limitando a auto-suficiencia da pro 
dugao agrTcola, geram o chamado mercado competitive entre as nagoes 
do mundo. Consequentemente, criam no panorama mundial dois grupos con 
correntes os paTses exportadores e os paTses importadores. 

Dentre os produtos basicos a al imentagao, o trigo tern des 
pertado grande interesse, tanto entre os paTses exportadores como entre 
os paTses importadores. No caso dos paTses exportadores como: Estados 
Unidos, Australia, Canada, Argentina, Franga, entre outros, o interesse 
esta relacionado com a colocagao desse cereal no mercado internacional 
atraves de pregos competitivos; enquanto os paTses importadores, tais 
como: Russia, China, Japao, Egito, Brasil, entre outros preocupam-se 
com a elevagao dos pregos do trigo no mercado mundial, procurando obter 
este cereal de modo mais economico para garantir o supriniento deste pro 
duto a populagao a pregos mais acessTveis. 

Urn exemplo do interesse pel a estimativa da produgao do 
trigoi por parte dos paTses exportadores, foi demons trado pel os Estados 
Unidos que empenharam grandes esforgos para estimar este cereal atraves 
de dados coletados por sateTite, ao iniciar em 1974 o projeto Large Area 
Crop Inventory Experiment (I.ACIE), cujo objetivo foi investigar a poten. 
cialidade do uso de sensoriamento remoto para urn monitoramento global 
da produgao agrTcola. 

0 cultivo do trigo no Brasil iniciou-se com os primeiros 
imigrantes europeus que se estabeleceram no Sul do PaTs. Atualmente, e£ 
te cereal e cultivado em todos os Estados da regiao Sul, Minas Gerais, 
Mato Grosso do Sul e Goias. Segundo Moura et alii (1979), a triticultu 
ra nacional vem experimentando sucessivas flutuagoes nestes ul timos anos 
em consequencia de adversidades climaticas como: geada em 1975, granizo, 
temperaturae umidade elevadas em 1976 e estiagem em 1977. 
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Mu 1 to embora suja meta do governo brasileiro tornar o 
PaTs auto-suficiente quanto a produgao do trigo, a demanda por cste ce 
real, no mercado interno, sempre foi maior do que a oferta, necessitate 
do desta forma, substanciais taxas de importagao. Atualmente, o trigo 
ocupa o segundo lugar na pauta de importagao do PaTs, superado somente 
pelo petrol eo. 


Diante desta situagao, pode-se concluir que a pressao 
que este cereal vem exercendo sobre o deficit do tesouro e da balanga 
comercial e bastante grande. Isto requer um previo conhecimento da nos. 
sa produgao, para uma tomada de decisao em tempo habil, garantindo as, 
sim o suprimento da populagao a pregos mais acessTveis. 

As informagoes que se tem sobre a area ocupada comacul_ 
tura do trigo vem sendo obtidas por metodos convencionais e, muitas v_e 
zes, sao subjetivas e insatisfatorias. 

Para fazer uma previsao de safras com confiabilidade 
prePfUf.alecida e em tempo habil, dois aspectos devem ser consider^ 
doai primeiro, fazer dentro da estagao de crescimento uma avaliagao p^ 
riodica, tanto da area como da condigao da cultura; segundo, uma previ 
sao da variagao mais provavel das condigoes de crescimento futuro e 
dos possTveis efeitos sobre a produgao. E extremamente difTcil obteres 
sas estimativas por metodos convencionais (MacDonald e Hall, 1978). Com 
relagao a avaliagao periodica, tanto da area quanto das condigoes da 
cultura, os dados coletados pelo sistema LANDSAT tern se mostrado uma 
ferramenta muito util devido as caracterTsticas multiespectrais e sind 
tica do sensor e de repetitividade do satellite. 

0 grandenumero de dados fornecidos pelo LANDSAT, dada 
a caractenstica de repetitividade, determinou a necessidade do desen 
volvimento de tecnicas computacionais para a interpretagao desses da^ 
dos visarido o monitoramento mais eficiente das culturas. Varios pes. 
quisadores como Bauer e Cipra (1973), Economy et alii (1974) e 
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Dietrich et alii (1975) afirmam que a analise dos dados do LANDSAT 
atraves de tecnicas computacionais e eficiente, rapida e de grande po 
tpncialidade. 


Dentro do contexto da analise digital supervisionada, urn 
grande problema e a utilizagao das estatTsticas de treinamento obtidas 
de pequenas areas para identificar e avaliar grandes areas agrTcolas. 
Como a culutra do trigo ocupa grande extensao de nosso territorio, mui^ 
tas variaveis como: epoca de plantio, tipos de solo, concentragao da 
cultura, si sterna de manejo, etc devem ser consideradas quando se pro. 
poe a utilizagao de estatTsticas de treinamento adquiridas em peque 
nas areas de trigo, para identificar e estimar outras. Neste caso, er 
ros significantes poderiam ocorrer devido a nao-representatividade das 
respostas espectrais da cultura do trigo, obtidas a partir dessas p£ 
quenas areas. Por outro lado, a cobertura aerofotografica de toda r£ 
giao tritTcola, aliada a informagoes de campo com a final idade de o_b 
ter todos os padroes de respostas espectrais para o treinamento de um 
Sistema Interativo de Analise Mul tiespectral , e posterior aval iagao dos 
resultados da classificagao, mostrar-se- ia impraticavel devido ao ele 
vado custo e tempo requerido para a analise. 

A amostragem tern demonstrado ser uma tecnica basica de 
coletar dados estatTsticos de agri cultura para um rapido e compreens_i_ 
vo sistema de informagoes exatas, atraves da tomada de amostras repre_ 
sentativas de somente uma pequena parte da area total. 

Em 1979 foram iniciados, no Instituto de Pesquisas Espa_ 
ciais (INPE), os primeiros estudos sobre a cultura do trigo, com obje 
tivo de verificar a potencial idade do uso de dados do LANDSAT para 
identificar e avaliar a area ocupada com esta cultura. De acordo com 
Mendonga et alii (1981), os dados coletados pelo LANDSAT mostraram-se 
potencialmente viaveis para a identificagao e aval iagao da area ocup£ 
da com esta cultura para as condigoes brasileiras. Dando conti nui'da'de 
a esta linha de pesquisa, o presente trabalho tern como objetivo princi^ 
pal desenvolver uma metodologia baseada em sistema de amostragem uti 1 i_ 
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zando dados de fotografias aereas e dados digitais do MSS do LANDSAT, 
para estimar a area ocupada com a cultura do trlgo. 

Para atender o objetivo proposto nesta pesquisa.varios 
aspectos foram analisados, tais como: selegao do classificador, sele 
gao da rnelhor maneira de obtengao das estatTsticas de treinamento pa_ 
ra o computador, selegao de um metodo de pos-processamento, determi_ 
rtagao do tamanho da unidade amostral e do tamanho otimo da amostra e, 
comparagao do desempenho obtido pelo sistema amostral com o metodo de 
cxpansao direta (somente dados de aeronave). Pel a impossibilidade de 
obter tais informagoes sobre toda area tritTcola brasileira, propo 
-se desenvolver este trabalho numa area pi 1 oto de 720 km , a qual ja 
dispunha de uma cobertura total de fotografias aereas. 



CAPITULO 2 


REVISAO BIBLIOGRAFIC A 

2.1 - anAlise de dados do landsat aplicados A agricultura 

As primeiras fotografias multi espectrais coletadas de cam 
pos agrTcolas foram adquiridas em 1964. Neste epoca ficou tambem reco 
nhecido o potencial dessas fotografias para identificar areas agrTcolas. 
Posteriormente, em 1967, foram implementados em computadores digitais me 
todos de reconhecimento de padroes, e uma das primeiras investigagoes 
usando Imagens adquiridas por satelite para identificar culturas foi 
executada por Anuta e MacDonald, 1974 (Bauer et alii, 1977). 

As tecnicas de sensoriamento remoto baseadas no uso de 
imagens orbitais, testadas em regioes agrTcolas, oferecem uma possTvel 
alternativa para substituir as fotografias convencionais, pois, alem de 
permitir obter informagoes do uso da terra mais frequentemente, devido 
a caracterTstica da repetitividade do satelite, podem tambem obter, in_ 
formagoes a cerca de areas mais extensas (Hay, 1974). 

Estudando quantitativamente a potencial idade dos dados 
MSS do LANDSAT e da analise digital para identificar e avaliar areas 
agrTcolas em tres municTpios ao norte do Illinois, Bauer et alii (1979) 
concluTram que os resultados obtidos com dados LANDSAT foram positivos, 
considerando o numero limitado de bandas espectrais em que atua o sensor 
MSS do LANDSAT e a baixa resol ugao espacial, quando comparados com os 
obtidos com dados de urn "Scanner" colocado a bordo de uma aeronave. Apia 
sar da baixa resol ugao espacial e do numero limitado de bandas espec 
trais, os dados adquiridos pelo LANDSAT nao apresentam problemas com o 
angulo de elevagao do sol, o que tern si do constatado em dados adquiri_ 
dos pelo "Scanner" aerotransportado. Continuando, os autores afirmamque 
outro fator que torna os resultados significativos e o tamanho da area 
levantada. 
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Para Hixson et alii (1981), a visao sinotica da terra 
fornecida por satelite de observagao, junto com o processamento digi 
tal dos dados, e muito util para identificar e estimar areas agncolas. 

Metodos de analise de dados obtidos pelo LANDSAT podem 
ser divi didos em: "interpretagao visual de imagens" e "tratamento nume 
rico ou digital dos dados gravados em fitas compatTveis com computador 
(CCTs)", No presente traba'ilio a analise dos dados do LANDSAT foi feita 
utilizando as tecnicas de interpretagao digital, atraves do 1-100. 
Por essa razao, o enfoque nesta segao foi mais direcionado a analise 
digital . 


Segundo Bauer (1976), na analise digital de dados de sen. 
soriamento remote, utilizam-se fundamental mente as variagoes espej: 
trais. Pode-se tambeni utilizar, se devidamente imp! ementadas, informa 
goes temporais e espaciais ern complementagao as informagoes provenien 
tes das variagoes espectrais. 

A obtengao de informagoes sobre as culturas, de acordo 
com Johnson e Coleman (1973), podera ser feita combi nando-se de forma 
integrada a capacidade do homem com a agilidade da maquina para obten. 
gao de resultados rapidos. 

Com relagao ao processamento digital, Hoffer ( 1 973) 
afirma que existem duas condigoes basicas para a definigao de uma clas_ 
se visando a analise de padroes dos dados de sensoriamento remote: a 
classe deve ser espectral mente separavel de todas as outras classes e 
de interesse do usuario. 

Para Higton e Bormann (1977), os erros na classificagao 
sao inevitaveis dada a complexidade da natureza. Eles afirmam que "e 
impossTvel identificar todos os tipos de alvos com base somente em suas 
refl ecfancias caracterTsticas na parte do visTvel e do infravemelho 
proximo do espectro eletromagnetico, sem associagao de criterios adj_ 
cionais tais como: tamanho, textura, forma e padrao". 
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Com relagao a confiabilidade dos dados do LANDSAT trata_ 
dos digitalmente, Economy et alii (1974) comentam quo o desempenho da 
classificagao nao depende somente da metodologia utiHzada parao reconhe 
clmento da resposta espectral, mas tambem da quail dade da informagao 
da verdade terrestre utilizada pelo analista para comparagao, 

Wigton e Steen (1973) visando a identificagao e avalia 
gao de areas ocupadas com culturas, a partir de dados do LANDSAT em 
areas selecionadas do Missouri, Kansas, Idaho 'e Dakota do Sul, explj_ 
cam o baixo desempenho da classificagao obtida (58,8^) como consequejn 
cia de grandes variagoes espectrais causadas pelo plantio das cultjj 
ras em diferentes solos. 

De acordo com Hlavka et alii (1980), urn problema funda_ 
mental na identificagao de areas agrTcolas com dados do LANDSAT e o nu 
mero de variaveis qua influenciam a reflectancia de uma determi nada 
cultura. Entre el as, tem-se o grau de maturidade da cultura que pode 
variar dentro de uma pequena area num dado tempo. Para os autores urn 
determinado campo pode estar no estagio de maturidade que estava o seu 
campo vizinho ha 18 dias atras. Nalepka et alii (1977), citados por e_s 
ses autores, afirmam que as diferengas observadas no estagio de desejn 
yolvimento do trigo no Estado de Kansas ocorriam nao so entre campos 
como tambem dentro de urn mesmo campo tritTcola. 

Ray et alii (1976) em pesquisas com dados de sensorU 
mento remoto, conduzidas em colaboragao com o pessoal do Centro Avanga_ 
do de Computagao (EAC) da Universidade de Illinois eo "Statistical 
porting Service (SRS)" do Departamento de Agricultura dos Estados Un_i_ 
dos (USDA), usando estimativa pel a regressao para estimar areas de cuj_ 
turas, concluiram que se os dados do LANDSAT gravados em fitas CCTs fo 
rem obtidos num tempo muito proximo a obtengao dos dados terrestres,i£ 
to e, sem defasagem de tempo, a estimativa da area agrTcola pode ser 
significativamente melhorada. 
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Utilizando os dados do LANDSAT-1, para analisar culturas 
agrTcolas e florestais em Michigan atraves do tratamento digital, Safir 
ot alii '(1973) concluiram que o reconheciniento de classes agrTcolas e 
florestais e bom desde que as copas ou folhagens das plantas apresentem 
um padrao uniforme, o que pode ser notado quando as plantas atingem a 
faso de maturagao. 

Confornie relata Hixson et alii (1980)., varios esquemas 
de analises digitais tern sido desenvol vidos e implementados em computa_ 
dores de varias Uni vers i dades , comercio e orgaos governamentais nos Es_ 
tados Unidos e em outros paTses. Cabe aos analistas de dados de senso_ 
riamento remoto determinar qual analise e mais apropriada ou qua! o al_ 
goritmo que melhor desempenha uma dada classificagao. 

Alem dos fatores como:diferentes variedades da cultura, 
tipo de solo, epoca de plantio, etc. que influenciam a resposta espe£ 
tral das culturas e, consequentemente, afetam os resul tados da classify 
cagao, para Hoffer (1977) outro fator que deve ser considerado e o trej_ 
namento do al goritmo no computador para reconhecer uma combi nagao parti_ 
cular de numeros representatives da reflectancia, em cada comprimento 
de onda, de urn alvo ou tipo de cobertura de interesse. 0 autor comenta 
ainda que uma classificagao digital de dados mul tiespectrais depende da 
definigao de urn conjunto de estatTsticas de treinamento, as quais sao 
na verdade representativas das caracterTsticas espectrais das varias 
feigoes da superfTcie da terra presentes nos dados da cena. 

Parsons e Fitzpatrick (1978) atraves de estudos para d£ 
terminar a melhor epoca de identificagao do trigo, em tres areas do di£ 
trito de Narrabri em New South Wales, Australia, atraves de interpret^ 
gao visual de imagens do LANDSAT-2, obtidas de maio a novembro, corre£ 
pondendo desde as epocas do plantio a colheita do trigo, concluiram que 
as epocas mais adequadas foram: junho, perTodo imediatamente apos a emer 
gencia, quando todas as areas plantadas aparecem como solo preparado 
devido a baixa densidade de cobertura do solo, e final de setembro,qua£ 
do a cultura do trigo esta formando a panTcula, porem,nao-atingindo ai£ 
da o-estagio de maturagao. 
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Para Hall et alii (1974), um grande problema durante a 
classificagao digital de culturas, usando dados de satelite, e a esco 
lha correta da area de treinamento para o computador, principalmente 
quando estas areas estao local izadas proximas ou na intersegao de dois 
tenias agrTcolas distintos mas de padroes espectrais bastante semelhari 
tes, Nesta situagao, recomendam os autores, "e conveniente nao tomar 
estas areas para treinamento". 

Com relagao a area de treinamento para a aquisigao de 
estatTsticas, Everitt et alii (1979) estudando a classificagao digi_ 
tal para fins de inventarios das terras ao sul do Texas, numa area de 
81.000 lia local izada nos municTpios de Kennedy e Willacy, utilizaram 
como treinamento do algoritmo no computador 0,5% da area total. 

Utilizando dados do LANDSAT-1 sobre tres areas ao norte 
do Illinois, Bauer e Cipra (1973) citam que as estatTsticas de treina_ 
mento, desenvol vidas para uma regiao, poderao ser usadas para classy 
ficar regioes adjacentes sem decrescimo do desempenho da classificagao 
digital. 


Bauer et alii (1979) obtiveram uma boa classificagao 
utilizando estatTsticas de treinamento, obtidas de areas distanciadas 
ate 80 km da area de estudo. Segundo esses autores, a distancia sobre 
a qual as estatTsticas de treinamento poderao ser estendidas com sjj 
cesso' dependera de uma serie de fatores como epoca do plantio, releva 
tipos de solo, etc.; porem, um importante requerimento e a similarida_ 
de entre a cena a ser classificada e a cena onde foram obtidas as est£ 
tTsticas de treinamento. 

Para Hixson et alii (1980), o procedi mento de treina_ 
mento e mais importante que a escolha do classificador para obter uma 
classificagao satisfatoria da cultura. Esses autores demonstraram em 
seus estudos que, fixando as areas de treinamento, houve pequena dif£ 
renga na exatidao da classificagao obtida quando se util iza ram diferen_ 
tes classificadores. 
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A definigao dos parametros de cada classe pode ser obti_ 
da atraves de abordagem supervisionada ou nao-supervisionada. Entretaji 
to, Wharton e Turney (1981) afirmam que em muitos casos a combi nagao 
das duas abordagens parece ser muito indicada, Esta combinagao, conti_ 
nuam os autores, e pratica para areas nas quals e difTcil ou impossTvel 
definir-se areas de treinamento para todas as feigoes, especificamente 
para aquelas que possuem areas pequenas, irregulares ou esparsamente 
distribuTdas. 


As experiences do Laboratory for Application of 
Agricultural Remote Sensing (LARS), de acordo com Hoffer (1973), tarn 
bem tem comprovado que a combi nagao destes dois sistemas (supervisio. 
nado e nao-supervisionado) mostram resultados mais satisfatorios. 

Depois de obtidos os parametros de cada classe presente 
» 

na area de interesse, seja pelo metodo supervisonado ou pelo nao-supe_r 
visionado, o restante da area e classificada utilizando o algorTtimo 
mais apropriado para a sua classificagao. Entretanto, os estudos feitos 
na area de sensor iamento remoto para identificar e avaliar areas agrT 
colas, utilizando a classificagao digital, tern demonstrado que os re 
sultados obtidos estao sempre sujeitos a erros de classificagao causa, 
dos pela semelhanga de respostas espectrais de alvos distintos. 

MacDonald e Hall (1978) comentain que a maior confusao 
entre as respostas espectrais das culturas e a resposta espectral do 
trigo ocorreu, principalmente, devido as culturas de cevada e centeio. 
Em regions onde essa confusao e acentuada, os trabalhoS'do LACIE d_e 
monstram que a estimativa da area de trigo.e diminuTcia devido ao acre£ 
cimo da area estimada com essas culturas. 

Hatfield (1981) analisando o comportamento espectral de 
82 variedades de trigo concluiu que houve pouca variagao espectral e_n 
tre as variedades ate a fase de emborrachamento, mas que esta variagao 
aumentou durante as fases de inchamento e maturagao. 
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Depois de feita a classificagao dos dados da area de i.n 
teresse, obtidos pelo LANDSAT, e sempre recomendavel que se faga a ava_ 
liagao de seu desempenho. Para avalia-lo, Hoffer (1973) comenta que as 
expert enci as do LARS tem demonstrado que uma combi nagao de tres tecni. 
cas: a) comparagao visual entre o "print-out" contendo os resultados 
da classificagao e o mapa tematico obtido a parti r da fotointerpreta. 
gao, b) amostragem de areas individual's designadas como areas testes, e 
c) analise da correlagao entre os dados adquiridos pelo LANDSAT e os 
dados de verdadc terrestre, fornecem a melhor indicagao sobre a exati_ 
dao da classificagao, 

De acordo com Chen et alii (1981), na ava liagao do d£ 
sempenho da classificagao devem ser considerados, conjuntamente, a pe_r 
centagem de classificagao correta e o erro de inclusao. 

2.2 - USD DO SISTEMA DE AMOSTRAGEM EM SENSORIAMENTO REMOTO PARA ESTIMAR 

Areas agrtcolas 


No Brasil, urn dos primeiros trabalhos utilizando fot_o 
grafias aereas como sistema de referencia a amostragem para estimar 
areas agrTcolas foi realizado por Schattan (1968). Neste estudo dados 
de fotografias aereas do MunicTpio de Tiete foram utilizados nos meto 
dos de expansao direta e de razao, para estimar areas cultivadas, pa_s 
tagem,,mata e outras classes de uso da terra daquele municTpio. 

Hay (1974) utilizando o metodo de regressao para estimar 
areas agrT colas em areas testes dos municTpios de San Joaquim (Califo_r 
nia) e Maricopa (Arizona) comenta que imagens orbttais e fotografias 
aereas utilizadas em conjunto podem aumentar a eficiencia e a exatidao 
da classificagao. 

Para Wigton e Bormann (1977), o uso de amostragem e 
mais economico para a estimativa de areas de umadeterminada cultura, 
principalmente em grandes areas onde e difTcil ou econoni camente i_n 
viavel urn conhecimento completo da regiao. Para estes autores, a exa_ 
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tidao da estimativa ira depender de quanto as caracterTsticas das amo£ 
tras sao representatives das caracterTsticas da populagao. 

Thomas e Hay (1977) em estudos cuja final idade foi l£ 
car amostras atraves de urn sistema de dupla amostragem, baseado nos dat 
dos do LANDSAT e nos dados do Statistic Report Service (SRS), do Depar 
tamento AgrTcola dos Estados Unidos (USDA), para serem util izadas na e£ 
timativa da area do trigo pelo metodo da regressao, concluiram que, se 
a correlagao entre as areas estimadas uti 1 1 zando dados do LANDSAT atr£ 
ves da analise digital e as mesmas areas obtidas por verdade terrestre 
for suficientemente grande, havera redugao na variancia da estimativa 
da area. 


De acordo com MacDonald e Hall (1978), a mais compree£ 
siva investigagao do uso de dados do MSS do LANDSAT, atraves de si$t£ 
ma de amostragem para levantamento de areas agrTcolas, principalmente, 
sobre a cultura do trigo, foi o Large Area Crop Inventory Experiment 
( LAC IE). Os objetivos do LACIE, comentam Hixson et alii (1981), foram 
assimilar a tecnologia atual de sensoriamento remoto dentro de urn si£ 
tema experimental e avaliar o seu potencial para determinar a produgao 
dentro de uma confiabil idade conhecida (limite de confianga de 90!?) e 
uma probabil idade de 1 0% de erro relativo para uma grande regiao ge£ 
grafica. A escala global do LACIE incorporou cerca de 15000 conjuhtos 
de dados adquiridos pelo LANDSAT-2, e mais de 2600 segmentos amostrais 
em 5 maiores regioes produtoras agrTcolas do mundo. 

Hanusckak et alii (1979), estimando areas cultivadas de 
mi 1 ho e soja em Iowa e trigo em Kansas, usaram como estimador das areas 
destas culturas dois metodos estatTsticos: expansao direta (somente d£ 
dos terrestres) e metodos de regressao (dados do LANDSAT e terrestres). 
Segundo esses autores, o ganho na precisao da estimativa da area obt£ 
da pelo uso da regressao quando comparada a expansao direta foi sub£ 
tancialmente maior; alem disto, este ganho tambem foi cons tatado nos e£ 
tudos feitos por Graig et alii (1979) ao estimar a area de trigo no E£ 
tado de Kansas . 



2.3 - A CULTURA DO TRIGO (Triviaum aootivwn L). 


0 trigo e uma planta herbacea,anual , monocotiledonea, da 
famTlia Graminoae, e do genero Triticum, 

Consta quo o Iraque foi o primeiro paTs no globo terres 
tre a conhecer o trigo ha cerca de 6500 anos atras. Porem, o primeiro 
paTs a cultiva-lo foi o Egito ha 5000 anos. Na China o trigo e cultiva. 
do desde 3000 anos atras. 

0 cultivo do trigo no Brasil ,segundo Moura et alii (1979), 
iniciou-se com os primeiros imigrantes europeus que se estabeleceram no 
Sul do PaTs. Segundo Almeida e Carnechio Filho (1973), a primeira notT 
cia que se tern da introdugao do trigo no Brasil data de 1556, quando 
se tentou o cultivo desta gramTnea na capitania de Sao Vicente. 0 Rio 
Grande do Sul destacou-se como principal produtor, atingindo em 1831 
expressiva area plantada. 

1) Estados brasileiros onde o trigo e mais cultivado. 

Rio Grande do Sul, como principal produtor, Parana, Sao Paulo, 
Santa Catarina, Mato Grosso do Sul e Minas Gerais. 

2) Clima 


A cultura do trigo, em escala mundial, tern a sua maior 
expressao em areas de clima temperado, do tipo mediterraneo, com inve_r 
nos umidos e veroes secos como os da Europa, centro dos Estados Unidos, 
Argentina, etc. 

De maneira geral , o trigo vegeta e frutifica bem nas re^ 
gioes em que, durante o ciclo vegetativo, a temperatura media e infe 
rior a 20°C e a precipitagao pluviometrica e superior a 180 mm, com boa 
distribuigao. 
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Com o melhoramento genetico, conta-se atualmente com cul^ 
tivarcs adaptadas aos climas tropicais, com veroes quentcs e umldos e 
invernos amenos e secos, As cultivares seleciwiadas para as condigoes 
brasileiras sao, cm goral * pouco exigentes quanto ao frio, indiferen 
tes ao fotoperiodismo, de ciclo curto, de boa arquitetura foliar» de 
grande tolerancia a umidade e de elevada resistencia as enfermidades. 

3) Solo 


0 trigo requer solos ferteis ,levemente acidos {pH entre 
5,5 e 6,5), com textura media e bem drenados. 

Segundo Almeida e Cornechio Filho (1973), de maneira g£ 
ral , o trigo nao pode ser semeado em terreno nem muito argiloso e nem 
excessivamente silicoso. Devem ser solos prof undos, para permitir a p£ 
netragao de suas raTzes que sao longas e descem verticalmente, e com 
boa capacidade de retengao de umidade. 

4) Variedndes 

0 genero Tvitiom compreende varias especies, distintas 
uma das outras, estando representadas por um grande numero de varieda. 
des, as quais se multiolicam por selegoes e cruzamentos. 

0 trigo, cultivado comercialmente no Brasil, pertence a 
especie Tritiaum aeotivum e a subespecie vulgara. Essa subespecie e 
constituTda por inumeras cultivares que estao em uso pelos agricultat 
res, outras em experimentagao e um numero maior de linhagens ou cultji. 
vares que foram superadas atraves do tempo e sao encontradas apenas em 
bancosde germoplasma como fonte de variabil idade genetica, para serem 
usadas em cruzamentos. 
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No caso do Brasil que possui uma ampla variagao edafocli. 
matica, ha necessidade do uma maior diversificagao das cultivares com 
caracterTsticas especTficas para cada regiao. 

Em geral, uma variedade de trigo para ser considerada 
boa precisa ter pelo menos tres caracterTsticas: ser produtiva.oferecer 
produto de boa qualidade e ser resistente as principais molestia. Alem 
disto, a variedade deve apresentar bons caracteres agronomi cos como 
adaptagao ao meio, resistencia ao degrane e ao acamamento. 

5) Epoca de plantio. 

A epoca de pi anti o esta condi cionada a tres fatores prin_ 
cipais; precipitagao pluviometrica, temperatura e epoca de colheita. 

0 trigo necessita de 60 a 90 di as de adequada umidade no 
solo durante o perTodo que vai da emergencia ao inTcio da maturagao, 

De urn modo geral, as cultivares existentes no Brasil sao 
resistentes a baixas temperaturas, sendo prejudicadas somente por ge£ 
das na epoca de floragao, Por outro lado, altas temperaturas em asso 
ciagao com dias longos reduzem significativamente o ciclo da planta e 
diminuem o potencial de produgao. 

A colheita deve ser feita em perTodo seco para a obte^n 
gao de graos de boa qualidade. 

Com base nestes fatores, deve-se fazer urn estudo para ca_ 
da regiao tritTcola, a fim de determinar a melhor epoca do plantio. 

6) Dcpnga e Pragas. 

Entre as principais doengas que causam serios prejuTzos 
economicos na produgao do trigo destacam-se: Ferrugem do col mo ( Puooi^ 

nia graminis tvitici ), Ferrugem das folhas ( Puooinia mbigovera 
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tioi) , helniintosporiose do trigo (HoUdnthoopovium oaHvum) f Mildlo ou 
Oidio do trigo (Eryaipho gvcuninio), 

Como principais pragas tem-se; lagarta mi 11 tar (LapJtygm 
frugipcrda) j lagarta do trigo (Pooudalctia adultora) i lagarta elasmo 
(Elaomopalpuo Hgnooollm)Q pul goes do trigo da ordom Homoptera e fa 
mTlia Aphididae. 

7) Colheita. 

A colheita do trigo e mecanizada, Depois de col hi do o 
trigo, e feita a secagem dos graos, em geral, com teor de umidade abaj^ 
xo de ]]%. Apos a secagem, o trigo e armazenado para a sua comercial^ 
zagao. 



CAPITULO 3 


MATERIAIS E METODOS 

3.! - AREA DE ESTUD O 

A cultura do trigo no Estado do Rio Grande do Sul encoii 
tra-se concentrada a noroeste do Estado e esparsamente distribuida ao 
sul . Para esse estudo utilizou-se uma area com aproximadamente 720 km 2 , 
na regiao de Cruz Alta, local izada dentro da area de concentragao de 

trigo, <entre os paralelos 28°30 ‘ e 28°45' S e meridianos 53 °30 1 e 

54°00 1 W (Figura 3.1) 


No municipio de Cruz Alta, de acordo com Brasil (1978), 
predominan as terras com uso mais apropriado para lavouras de ciclo 
curto ou longo, permitindo em quase sua totalidade dois cultivos por 
ano, quanto a fertilidade apresentam serios problemas, pois 70% de 
suas terras necessita, de altas aplicagoes de fertilizantes e correti_ 
vos de acidez do solo. Esta regiao e tambem caracterizada pel a grande 
exigencia de pratica conservacionistas; porem, trata-se de grandes e_x 
tensoes de terras proprias a mecanizagao. 

Basicamente tres fatores contribuiram diretamente para 
a escolha da area de estudo no municipio de Cruz Alta: a) a grande coji 
centragao da cultura do trigo; b) a representatiVidade da cultura a nT 
vel estadual e; c) a existencia de informagoes preliminares oferecidas 
por entidades governamentais e privadas que atuam na regiao. 

0 municTpio de Cruz Alta possui uma area de 3.462 km 2 
e aareade estudo corresponde a cerca de 23% deste municipio. 
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ORIGINAL PAGE 13 
OF POOR QUALITY 



Fig. 3.1 - fire a de concentragao da cultura do trigo no Estado do Rio 
Grande do Sul . 


FONTE: Brasil (1978). 
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3.1.1 - CLIMA 


Mota, citado em Brasil (1973), classificou o clima do Es 
tado do Rio Grande do Sul no sistema de Koppen e verificou a ocorren 
cl a de duas variedades cl 1 mat leas especTficas: 

a) "Cfa" clima subtropical (ou virglnlano) umido sem estlagem .A 

teniperatura media do mes mais quente e superior a 22 °C 
e a do mes mais frio varia de 3 a 18°C. 

b) "Cfb" clima temperado em que a temperatura media do mes mais 

quente e inferior a 22°C. 

0 municTpio de Cruz Alta, segundo Mota, encontra-se sob 
o regime de clima temperado (Cfb) com lima temperatura media anual de 

180c. 

3.1.2 - PRECIPITAQAQ 


No Rio Grande do Sul, a precipitagao media anual varia 
usualmente de 1100 a 2400 mm e esta relacionada as diferentes regioes 
f isiograficas . As chuvas sao muito variaveis e esta variagao e verify 
cada nao so de urn ano para outro, mas tambem durante urn mesmo ano. Na 
maioria dos anos, os meses mais chuvosos sao*, maio, junho e setembro e 
os menos chuvosos novembro, dezembro.e fevereiro, (Brasil, 1973). 0 mu 
nicipio de Cruz Alta situa-se numa regiao onde a precipitagao pluviom£ 
trica esta entre 1750 e 2.000 mm. 

3.1.3 -SOLOS 


Segundo o Levantamento de Reconhecimento dos Solos do 
Estado do Rio Grande do Sul, no municTpio de Cruz Alta ocorrem solos 
de tres Unidades de mapeamento: Cruz Alta, Santo Angelo e Tupancireta 
(Brasil, 1978). 
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1 ) Unidade de Mapeamento Cruz Alta 


E constituTda por Latossolo Vermel ho Escuro Distrofico, 
textura media e substrato arenito. A area abrangida por esta unidade 
em todo o Estado perfaz um total aproximado de 7945 km 2 , o que repre_ 
senta cerca de 3$ da area territorial do Rio Grande do Sul. Seu relevo 
geral e ondulado, constituTdo por elevagoes arredondadas, que forrram 
depressoes fechadas. 

2) Unidade de Mapeamento Santo Angelo 


E formada por Latossolo Roxo Distrofico, textura argilo_ 
sa, situada em relevo ondulado a suavemente ondulado, formado por ele 
vagoes com centenas de metros de comprimento, cuja declividade varia 
de 3 a 10$, podendo ocorrer relevo mais forte com, no maximo, 15?$. 

3) Unidade de Mapeamento Tupancireta 


E formada de Podzol ico Vermel ho Amarelo, textura media, 
ocupando area de aproxiniadamente 795 km 2 , o que representa 0,3$ da area 
do Estado. Esta unidade e encontrada nos municipios de Tupancireta, JO 
lio de Castilho e Cruz Alta. 

3.1.4 - VEGETAQflO 


Sobre a area de estudo encontra-se como vegetagao natu 
ral o campo, na maioria das vezes utilizado como pastagem, pequenas 
areas com matas que aparecem de forma isolada e mata de galeria que 
acompanha os cursos d'agua. Mas depressoes a pastagem e um pouco m£ 
lhor, composta por diversas gramTneas. 



3.2 - LEVANTAMENTO AEREO 
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Em 2 de setembro do 1979, realizou-se sobre a area de es. 
tudo um sobrevoo utilizando camera metrica do tipo RC-10, com distancia 
focal de 152 mm e films Aerochrome 2443, cuja sensibil idade espectral vai 
de 400 a 900 nm. Utilizou-se ainda urn filtro com a finalidade de bio 
quear toda a luz com comprimento de onda na fnixa do azul , diminuindo ajs 
sim o efeito atmosferico. 0 voo foi realizado com recobrimento lateral 
e longitudinal de 30$, perfazendo 7 faixas de voo num total del04fotos. 
Transparencies positivas foram obtidas na escala de 1:20.000. 

As fotografias aereas foram plotadas em cartas topografi_ 
cas do IBGE na escala de 1:50.000 para facilitar o trabalho de localize 
gao dos talhoes de trigo no campo. 

3.3 - AQUISIQAO DOS DADOS LANDSAT 


A epoca ideal para a aquisigao dos dados LANDSAT para a 
identificagao de trigo e durante o estagio de maturagao (Chen et alii, 
1981), porque nesta epoca o trigo apresenta uma coloragao amarelo-ouro di_ 
ferenciando-se das outras culturas (predominantemente pastagem) que ain_ 
da permanecem verdes. Para aquela regiao, no ano de 1979, a epoca ideal 
seria final de setembro ou inTcio de outubro. Entretanto, como os dados 
LANDSAT do dia 22 de setembro apresentaram 100% de cobertura de nuvens, 
tornou-se impossTvel a sua utilizagao. Para este estudo foram utilizados 
dados LANDSAT, gravados em fitas compatTveis com c;mputadores (CCTs), r£ 
ferentes a orbita 220, ponto 32 da passagem do dia 4 de setembro de 1979. 
Nesta epoca, a maior parte do trigo encontrava-se em estagio de emborra. 
chamento. 
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3.4 - SISTEMA UTIL IZAPO NO TRATAMENTO DIGITAL DOS DADOS DO LANDSAT 

0 Si sterna Interativo de Analise de Imacjeni Multiespe£ 
tral, IMAGE-100 ou simplemente 1-100, produzido pela General Eletric, 
tem side utilizado com a fungao principal de extrair informagdes tenia, 
ticas das imagens multiespectrais. 

0 si sterna permite urn treinamento intehativo pelo anal is. 
ta (classificagao supervisionada) ou de maneira sisteniatica, sem parti_ 
cipagao direta do analista (classificagao nao-supervisionada) , identify 
cando-se sobre a area de estudo pontos com caracterT sticas semelhantes. 
Ao conjunto de pontos com caractensticas semelhantes atribui-se uma 
cor, urn dTgito ou urn valor numerico, dependendo do local de saida dos 
resultados, que podera ser atraves de urn terminal grafico alfanumerico, 
um visual izador colorido e uma impressora de linha, 

Caractensticas e potencial idades do sistema 1-100 (Ge 
neral Eletric Company, 1975): 

a) permite a analise do comportamento espectral ate oito temas dji. 
ferentes simultaneamente; 

b) permite a analise das imagens em tempo real, possibil itanto a 
interagao direta anal ista/imagem que e fundamental para a o.b 
tengao de resultados de alto graU de confiabil idade; 

c) utiliza tecnicas discriminatorias com parametros facilmentem£ 
dificaveis pelo analista, em tempo real, para a eliminagao da 
classificagao indesejada; 

d) permite a aplicagao de fungoes de pre e pos-processamento so 
bre a imagem, cuja final idade e melhorar sua qual idade para ob. 
tengao das informagdes necessarias a analise digital; 
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e) permite varios tipos de saTdas dos resultados da classificagao 
como: visual izador colorido, terminal grafico alfanumerico, 1m 
pressora de linha e fitas magneticas para posterior obtengao 

* 

de resultados em imagem (ou mapas). 

3.5 - METODOLOGIA 

3.5.1 - INTERPRETACAO VISUAL DAS FOTOGRAFIAS AEREAS 

De posse das fotografias aereas, em transparent ias po. 
sitivas, procurou-se fazer urn estudo preliminar, auxiliado pelos co 
nhecimentos de campo, a fim de reconhecer as caracterTsticas tonais 
dos diferentes alvos de ocupagao do solo. Posteriormentejestabeleceij 
~se uma chave de interpretagao cuja final idade foi ajudar ofOtointer 
prete a orgahizar as informagoes presentes nas fotografias aereas. 

Durante a fotointerpretagao procurou-se caradterizar 
os alvos de ocupagao do solo, dividindo-os em subclasses de acordo' 
com os padroes de tonal idades apresentados nas fotos aereas. Esta di_ 
visao teve o objetivo de auxiliar o entendimento do comportamento es_ 
pectral destes alvos quando observados nos dados do MSS do LANDSAT. 

Alem do criterio basico da cor, utilizado para cara£ 
terizar os diferentes temas, usou-se tambem os atributos textura foto 
grafica e arranjo espacial para auxiliar as decisoes sobre interpre 
tagao. Assim, todas as feigoes contidas nas fotos puderam ser identi_ 
ficadas corretamente. 

Os resultados finais da fotointerpretagao serviram p_a 
ra verificar a eficiencia do metodo durante o desenvolvimento da pe_s 
quisa e como uma das variaveis utilizadas na estimativa da area tri_ 
tTcola. 
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3.5.2 - ANALISE DIGITAL DOS DADOS DO LANDSAT 


Para a analise digital do trigo, foram utilizadas info_r 
magoes dos 4 canais do MSS do LANDSAT, adquiridas em uma unica data e 
gravadas em fitas CCTs. 

A fim de manter uma relagao de 1:1 entre os elementos 
de resolugao dos dados do LANDSAT e do visualizador do I -1 00 , fez-se 
uma ampliagao da cena correspcndente a area de estudo para a escala de 
1 : 100 . 000 . 


Como a area de estudo tem uin formato retangular de apro 
ximadamente 40 km x 20 km {40 x 20 cm na escala de 1:100.000) e devi_ 
do ao visualizador do I -100 apresentar uma area quadrada de 30 x 
30 cm, a area foi dividida em dois modulos (A e B), cada urn com uma 
area aproximada de 20 x 20 km. Para localizar os vertices da area de e^ 
tudo, tomou-se, com o auxTlio do cursor, umaregiao que supostamente 
continha urn dos vertices de cada modulo, a qual foi ampliada para a 
escala de 1:20.000, compatTvel com a escala do mapa tematico obtido 
atraves da fotointepretagao. A comparagao visual das informagoes, tan_ 
to no mapa tematico quanto no visualizador do 1-100, permit! u ao ana 
lista julgar e decidir sobre cada vertice e, com o programa trago, fez^ 
-se a dolimitagao definitiva de cada modulo. 

Apos a delimitagao da area de estudo, fez-se a clas. 
stficagao digital utilizando o metodo HTbrido M-2, e como pos-processa_ 
mento fez-se uma filtragem especial utilizando o metodo de Uniformiza 
gao de Temas-opgao UNIT0T com limiar 2 e peso 2, descritos a seguir. 
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Metodo HTbrido (M-2) 

Este metodo utiliza, inicialmente, uma abordagem nao-s_u 
pervisionada (K-MtDIAS) para a separagao das feigoes da area de estudo 
em classes espectrais homogeneas, Posteriormente, areas de treinainento 
com "pixels" puros {"pixels" localizados nos centros dos talhoes sem in. 
fluencia de falhas de cultura) sao adquiridas destas classes espectrais 
para a obtengao das estatTsticas de treinainento utilizadas na abordagem 
supervisionada (MAXVER). Para cada abordagem sera descrito o classifica 
dor envoi vido , 

1) K-MEDIAS 

0 classificador K-MEDIAS define classes espectrais a 
partir de urn conjunto de celas em quatro dimensoes. 0 analista fornece 
ao sistema o ntimero de centros iniciais (vetor media inicial para cada 
classe), a precisao e o numero de iteragoes que deseja, 0 sistema divi_ 
de os dados da area de estudo associando cada vetor ao centro mais pr5 
ximo, segundo a fungao discriminante distancia euclidiana. Apos esta a_s 
sociagao, sao calculados novos centros, Este ciclo repete-se ate que: 

a) atinja o numero de iteragoes desejado, ou 

•b) a diferenga entre o centro anterior e o atual nao seja maior 
que a precisao escolhida. 

Extrai-se entao o vetor media e a matriz de covariancia 
para cada grupo de dados que serao, posteriormente, utilizados na cla£ 
sificagao. 


Em virtude do tamanho relativamente grande da area e da 
utilizagao das informagoes dos quatros canais do MSS do LANDSAT, houve 
problema de restrigao de memoria contputacional para armazenar inform^ 
goes de todos os "pixels" para o processo de agrupamento; desta forma, 
o algorTtmo de agrupamento K-MEDIAS foi utilizado para agrupar "pixels" 
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resultantes de uma amostragem sistematica, Tal amostragem foi '■ felta 
usando uma grade cujas llnhas e colunas foratn tragadas a cada 10 "pi 
xels" da imagem; todos os "pixels" da area de estudo que se situavam 
nos pontos de intersegao das linlias e colunas foram utilizados na amos 
tragem. 

2) SISTEMA DE CLASS IF I CAQAO GAUSSIANA (MAXVER) 

0 MAXVER e uni si sterna supervisionado que se destina a 
classificagao ponto a ponto de imagens multiespectrais. 0 criterio uti_ 
lizado na classificagao e o de maxima verossimilhanga, segundo clas_ 
ses escolhidas interativamente pelo usuario, 0 uso do sistema MAXVER 
divide-se em duas fases distintas: a) fase de treinamento, na qua! o 
usuario escolhe as classes e adquire.atraves do cursor, areas de trej_ 
namento para obtengao das estatTsticas de treinamento (vetor media e 
matriz de covariancia) e; 2) fase de classificagao propriamente dita, 
onde o restante da imagem e classificoda da seguinte forma: para urn 
dado ponto X qualquer, sao calculados os valores de probabil idade de_s 
te ponto pertencer a cada uma das classes definidas na fase de tre_i_ 
namento. 0 ponto X e entao associado a classe'cuja probabil idade e 
maior. 


0 sistema MAXVER tern a capacidade de analisar ate 32 
classes com urn numero de 160 posigoes de tomada de area de treinamen^ 
to, Descrigoes mais detalhadas deste sistema podem ser encontradas em 
Velasco et alii (1978), 

Uniformizagao de Temas-opgao UNITOT 

0 UNITOT e utn dos metodos de pos-processamento que cort 
siste em uma filtragem espacial, o qua! analisa o "pixel" central de 
uma matriz de 3 x 3 "pixels". Para a analise, doi s parametros devem 
ser considerados: 1 i mi ar e peso. Estes sao predeterminados pelo anal is. 
ta. 
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A frequencia de ocorrencia do cada classe no interior 
da niatriz e determinada. A classe do "pixel" central e contada taji 
tps vezes quanto for o valor do peso. Posteriormente, o valor de fre 
quencia da classe de maior ocorrencia no interior da niatriz e compa. 
rada com o limiar. Se o valor do limiar for menor, a identidade do 
"pixel" central sera substituTda pel a classe de maior ocorrencia deji 
tro da niatriz. Se o valor do limiar for maior, a classe do "pixel" 
central permanecera sem alteragao. Desta forma o algorTtmo compare 
todos os "pixels" contidos na area de estudo (Dutra, no prelo). 

0 uso dos dados digital izados do MSS do LANDSAT como 
uma variavel nos levantamentos agncolas, de acordo com Sigman et 
alii (1977), produz erros de estimativa que podem ser divididos eni 
duas categorias quanto as origens: a) erros procedentes do classify 
cador devido as diferentes fungoes descriminantes (linear ou quadra, 
tica), ao conjunto de treinamento, a probabi 1 idade a priori e ao nD 
mero de categorias de classificagao usada no delineamentg da classi_ 
ficagao; b) erros de niedidas devido ao tamanho do elemento de resold 
gao do LANDSAT ("pixels"^ 57 x 79 m). 

Quanto aos erros de classificagao, estes podem ser mi_ 
nimizados atraves da escolha do metodo apropriado durante a analise 
digital, ao passo que, com relagao a segunda fonte de erros, estes 
sao inerentes ao si sterna sensor utilizado e so podem ser minimize 
dos pelo avango tecnologico de futuros satelites com melhor resol u_ 
gao. 


A razao da utilizagao do metodo HTbrido M-2 e do pos- 
processamento UNITOT foi para procurar minimizar os erros inerentes 
a classificagao digital. Estes aspectos foram, anteriormente, abord^ 
dos em tres trabalhos sucintamente descritos a seguir. 
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Escolha do melhor classificador para identificar e avaliar a area 

I ' I l « l >■> l »l >■* I l * ■ l» n ■ I l* ■ fcw * « «l I ■ ■ i • » «' » * • • » • " » J » ' ■ 'll — — '■ 

de trigo. 


Este aspecto foi previamente estudado e descrito por Chen 
et alii (1981), que realizaram estudos para comparar o desempenho da 
classificagao do trigo utilizando os classificadores: metodo do parale 
lepTpedo (SINGLE-CELL), metodo de classificagao nao-supervi s i ona da 
("clustering") K-MEDI AS e metodo de classificagao gaussiana (MAXVER) . Pa^ 
ra os classificadores MAXVER e K-MEDI AS , foram utilizados dois metodos 
diferentes para adquirir as areas de treinamento. Quanto as areas de 
treinamento para o trigo foram divididas em quatro subclasses, adquirj^ 
das subjetivamente pelo analista de acordo com os diferentes padroes t£ 
nais observados (em falsa cor) no visual izador do 1-100, o classifies 
dor MAXVER foi o metodo que mostrou melhores resultados com alta perce^n 
tagem de classificagao correta (87,3%), erros de inclusao relati vamente 
baixo (1 2,9%) e menor diferenga relativa (+ 0,2%). 0 metodo do parale. 
lepTpedo tambem foi testado com opgao MULTI-CELL (General Eletric } 1975) 
porem os resultados foram muito inferiores em relagao aos apresentados 
na Tabela 3.1. 
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TAOELA 3.1 


DESEMPENHO DA CLASS IFICAQAO PONTO A PONTO PARA VflRIOS 
CLASSIFICADORES E PROCEDIMENTOS DE TREINAMENTO 


CLASSIFICADOR 

— 
DESEMPENHO DA CLASSIFICAgXO(%) 

CLASS IFICAgXO 
CORRETA CCC) 

ERRO DE IN 
CLUSA0 (El) 

DIFERENgA 

RELATIVA 

SINGLE-CELL 

88,1 

15,8 

-4,5 

MAXVER (a) 1 

84,6 

14,8 

-0,7 

MAXVER (b) 2 

87,3 

12,9 

+0,2 

K-MEDIAS (a) 3 4 

88,7 

16,0 

+5,6 

K-MEDIAS (b)- 

75,3 

10,7 

-15,7 


1 (a) As areas de treiriamento obtidaa das diferantea claaoea de tri 

' go formaram uma mica olaooe para a aquiaigao das eatatiatZ 
aae do treinamento (rotor media a rntriz de covarianaia) utZ 
lizadas durante a claaaifiaagdo . 

2 (b) Para eada classe de trigo 3 obteve-ae uma estatistica > de trei 

namento. 

3 (a) Aa claaoea form agrupadaa aegundo o criterio nao-interativo. 

4 (b) Vtilizouree pare, inicializagao do "clustering" apenaa as in 

formagoea contidaa naa areas de treinamento obtidas para os 
outroa claaaificadjorea. 
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Escolha do metodo de agulsigao das estatTstieas de treinamento. 

Uma vez definido o classificador , outro aspecto levado em 
consideragao foi a escolha do melhor processo de aquisigao das estatTs 
cas de treinamento (vetor media e matriz de covariancia). Este aspecto 
foi abordado por Lima et alii (1982) que compararam os rosultados da 
classificagao do trigo, obtidos utilizando dois metodos de classifies 
gao digital; o nietodo M-l (MAXVER-b), onde as estatTsticas de treinamen. 
to sao adquiridas em talhoes conhecidos, representati vos dos alvosdein 
teresse, e o nietodo Hibrido, o qual foi chamado metodo M-2 que consiste, 
inicialmente, na utilizagao do classificador nao-supervisionado K-ME 
DIAS ("clustering"), onde, atraves da fungao discriminante (Distancia 
Euclidiana), faz-se o grupaniento dos "pixels" de respostas espectrais sj_ 
ml T ares . Posteriormente, as areas de treinamento com "pixels" puros fo 
ram adquiridas e utilizadas pelo MAXVER. 

Esses dois metodos de aquisigao das estatisticas de trei_ 

namento foram aplicados sobre a mesma area de estudo e, posteriormente, 

o desempenho de cada metodo foi aval i ado com base em alguns fatores, 

tais como a classificagao correta (CC), o erro de inclusao (El) e a di_ 

ferenga relativa (DR), utilizando-se as areas estimadas atraves das fo 

tografias aereas e do 1-100. Para este estudo, utilizaram-se os result^ 

dos da classificagao digital e da fotointerpretagao de cinco segmentos 
2 

de 40 km cada, com diferentes densidades de culture do trigo. 

Para comparar o desempenho dos dois metodos, aplicou-se 
o teste t (a=0,05) para cada urn dos fatores mencionados no paragrafo an. 
terior, o que nao revelou diferenga significativa. 0 uso do metodo HT_ 
brido M-2 para obter as estatTsticas de treinamento, entretanto, requer 
menos interagao homem/maquina, o que torna o metodo menos subjetivo. 
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Filtragem espacial 


0 terceiro aspecto considerado foi a escolha de um meto 
do de pos-processamento, baseado numa flltragem espacial, para refinar 
os resultados da classificagao. Este aspecto foi abordado por Moreira 
et alii (1982), que utilizaram como pos-processamento o metodo de Uni_ 
formizagao de Temas-opgao UNITOT. Este metodo envoi ve dois parametros 
determinados pelo analista: limiar e peso. Neste estudo foram testadas 
as seguintes combi nagoes de limiares e pesos: (1,2); (2,1); (2,2); 
(2,3); (3,2); (3,4) e (3,5). Os pontos classificados conio trigoparaca. 
da combi nagao foram impressos em forma de mapa alfanumerico na escala 
de 1:20.000. Fez-se a comparagao de cada mapa sobrepondo-se a cada uni 
deles o mapa tematico obtido a parti r da fotointerpretagao. Esta com 
paragao foi feita visual mente, sobre uma mesa de luz, observando-se os 
erros de omissao e de inclusao. Com base nesta comparagao, a conibinagao 
limiar 2 e peso 2 foi considerada a melhor. 

Para analisar quantitativamente o efeito da combi nagao 
limiar 2 e peso 2 sobre os resultados da classificagao, os autores aci_ 
ma citados fizeram analises estatTsticas selecionando cinco areas tes_ 
tes (~ 40 km 2 cada) com diferentes densidades de trigo, tanto no mapa 
alfanumerico adquirido por esta combinagao, quanto no mapa tematico. Pa_ 
ra este estudo, aplicou-se o teste t com c t = 0,05 aos dados de: classi_ 
ficagao correta (CC), erros de inclusao (El) e areas estimadas. 0 r£ 
sultado do teste mostrou que o uso do UNITOT melhorou significativame_n 
te a classificagao correta e diminuiu o erro de inclusao. 

Posteriormente, com base nestas areas testes ,esses auto_ 
res calcularam a diferenga relativa e a raTz media quadratica dos -e^r 
ros. Com a filtragem espacial, a diferenga relativa mostrou uma redu_ 
gao de 4,5% para 2,7%, e a raTz media quadratica dos erros diminuiu de 
126,18 ha para 107,02 ha. Com base nestes resultados, os autores, con_ 
cluiram que a aplicagao do metodo UNITOT oferece uma estimati va da area 
mais precisa do que aquelas obtidas sem a sua utilizagao, Porem, pelo 
baixo erro de inclusao e maior classificagao correta, este metodo aiu 
mentou a exatidao de mapeamento, conforme pode se observar na Tabela 3.2. 
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TABELA 3,2 

DADOS DE CLASS I FI CACAO CORRETA E ERRO DE 
INCLUSftO PARA AS A R EAS TESTES 


AREAS 

TESTES 

CLASS IFICAQAO CORRETA 
(CCJ5) 

ERRO DE 
(El* 

INCLUSAO 

S/UNITOT 

C/UNITOT 

S/UNITOT 

C/UNITOT 


1 

83,75 

86,05 

14,73 

12,26 

2 

82,83 

83,93 

17,73 

14,67 

3 

77,51 

78,93 

22,89 

20,23 

4 

88,41 

89,08 

28,4b 

23,98 

5 

87,85 

88,30 

21,60 

21,15 


Coni base nestes estudos, definiu-se que seria uti 1 i zado 
neste trabalho o metodo Hibrido M-2 e como filtragem espacial, o meto 
do de Uniformizagao de Temas-opgao UNITOT, com limiar 2 e peso 2. 

3 . 5.3 “ estimativa da Area qcupada com trigo, pelq usd de amostragem, 

UTILIZANDO DADOS DO MSS DO LANDSAT 


Esta fase consistiu em avaliar a viabilidade do uso do 
processo de estimate va por amostragem para o levantamento da area de 
trigo. 


Dentre os metodos estatTsticos existentes, destinados 
aos levantamentos por amostragem, dois deles foram selecionados para 
o estudo: a) estimativa por Tndice e, b) estimativa pel a regressao. Am 
bos os metodos podem fornecer uma boa precisao da estimativa da area de 
trigo, se a variavel auxiliar X (i.e., area estimada atraves dos dados 
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do LANDSAT) for correlacionada com a variavel Y (i.e., estimativa de 
area obtida por fotografias aereas ou campo). Qaunto maior foracorre 
lagao entre as variaveis, menor a variancia da estimativa pela regres 
sao (Hanuschack et alii, 1979). 

Quando a relagao entre X e Y for aproximadamente linear 
e a reta que a representa passar pela origem dos eixos, oprocessodos 
indices e mais apropriado. Se a reta que representa esta relagao nao 
passar pela origem dos eixos, o processo de regressao e o mais i nd i c a_ 
do (Cochran, 1965). 

0 processo estatTstico empregado nesta fase foi dividi_ 
do em duas etapas. A primeira etapa consistiu no estudo da correlagao 
existente entre as areas de trigo, obtidas a partir dos dados do MSS 
do LANDSAT com as mesmas areas estimadas atraves de fotos aereas, va_ 
riando-se o tamanho dos segmentos. Com os resultados da anal ise decor 
relagao, foi possTvel verificar quais tamanhos de segmentos deveriam 
ser utilizados para testar os metodos de estimativa por amostragem. A 
area de estudo foi dividida em 72 unidades de aproximadamente 10 km 2 . 
Testaram-se os tamanhos de segmentos (unidades amostrais) de 10, 20, 

30, 40 e 60 km 2 . Para tamanhos de segmentos acima de 10 km 2 j as areas 
foram formadas pela aglutinagao de 2, 3, 4 e 6 unidades de 10 km 2 do 
seguinte modo: 

1) para os segmentos com areas de 20 e 30 km 2 , tomaram-se unida^ 
des adjacentes no sentido horizontal, vertical e diagonal da 
area de estudo; 

2) para segmentos com areas de 40 kni 2 , tomaram-se unidades adja_ 
centes de areas de 2 x 2 unidades no sentido horizontal e ver 
tical ; 
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3) para tamanho do segmentos de 60 km 2 , usou~$e 3x2 segmentos no 
sentido horizontal e 2 x 3 segmentos no sentido vertical. 

Devido a area de estude ter urn tamanho fixo de aproximada 
niente 720 km 2 e como foram testados diferentes tamanhos de segmentos, o 
numero total de unidades da populagao, para cada tamanho de segmentos, va 
riou con forme a Tabela 3.3. 


TABELA 3.3 

N0MERO de unidades amostrais da populace para os 

DIFERENTES TAMANHOS DE SEGMENTOS TESTADOS 


TAMANHO DO 

SEGMENTO (,km 2 l 

10 

20 

30 

40 

60 

N9 DE UNIDADES 
AMOSTRAIS DA 
POPULAQftO 

72 

36 

24 

18 

12 


As areas estimadas atraves da interpretagao visual das 
fotos aereas e da analise digital dos dados do MSS do LANDSAT, eni cada 

segmentos, foram obtidas da seguinte maneira: 

a) sobre o mapa alfanumerico obtido na etapa de classificagao dj_ 
gital, fez-se a contagem manual das areas: (1 ) classificadas cor 
retamente como trigo; (2) de nao-trigo classificadas como trigo 
e; (3) de trigo, nao-classificadas ; 

b) . as areas ocupadas com a cultura do trigo, estimadas atraves do 

si sterna 1-1 00 utilizando dados do MSS do LANDSAT, foram deterim 
nadas somando as areas classificadas conforme os subitens 1 e 
2 do item (a) ; 
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c) as areas de trigo obtidas atraves da interpretagao visual de f<a 
tografias aereas foram determinadas somando as areas classifies 
. das segundo os subitens 1 e 3 do item (a). 

Desta forma, varias tabelas contendo as areas desta cul_ 
tura, obtidas atraves de fotografias aereas e analise digital, foram g£ 
radas de acordo com o tamanho do segmento e o sentido (horizontal, ver 
tical e diagonal) considerado na area de estudo. Com os dados contidos 
nestas tabelas fez-se uma analise de correlagao. 

Inicialmente estes dados foram plotados num si sterna de 
coordenadas cartesianas, obtendo-se para cada tamanho de segmento uni 

diagrama de dispersao, Obtidos os diagramas de dispersao, notou-se uma 
relagao linear entre as variaveis. Estes dados foram submetidos a uma 
analise de regressao linear simples, onde foi empregado o metodo dos nn_ 
nimos quadrados. Apos o estabelecimento da regressao linear simples, pr£ 
curou-se verificar se os coeficientes lineares das retas diferiram e£ 
tatTcamente de zero (a* 0,05). 

Tendo em vista que a reta obtida para cada conjunto de 
dados nao passou pel a origeni e que os dados do satelite e aeronave apre_ 
sentaram uma alta correlagao, conforme pode ser vista na Tabela 4.1, o£ 
tou-se pelo uso do metodo de estimativa pel a regressao, cuja formula m£ 
temat'ica e descrita a seguir (Cochran, 1965), 

A estimativa da area de trigo utilizando esta abordagem 

e obtida por: 

? R = N. Y R , 

onde: 


N = numero de segmentos da populagao. 
? = area total estimada com trigo. 
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Y r = y + b (X - x). 
n 

y « E y./n - area media de trigo obtida dos dados de fotografias 
aereas para (n) unidades amostrada. 


5< = E x^/N - area media de trigo obtida atraves dos dados do 
LANDSAT para toda area de estudo {N unidades). 

n 

x « E x^/n - area media de trigo obtida dos dados LANDSAT para 
(n) unidades amostradas. 

. f (x, - X ) . (y i - y) 
b 

f (x, - x) 2 


onde: x^, y . , etc. sao definidos como aci 


ma. 


A estimativa da variancia para pequenas amostras e dada 


por: 


V(T > . ■£ . JLiJl . S(y. -y) 2 .L^ 2 . { , + l +n (Lzil 2 , 

" " ‘ » - 2 11 t (Xj - X ) 2 


onde: 


n 


{f (y i - y) . (x i - x) } 2 
( y i - y) 2 . (x i - X) 2 } 


r 2 - coeficiente de determinagao entre as areas de trigo obtidas 
atraves das fotografias aereas e dos dados do LANDSAT para as 
amostras selecionadas. 
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Na segunda etapa utilizaram-se, na equagao acima, os 
dados obtidos nas fases anteriores, de interpretagao visual das fotos 
aqreas e da analise digital dos dados do LANDSAT, para estimar a area 
ocupada coni trigo. Para que o emprego deste sistema de amostragem possa 
fornecer uma estimativa da area de trigo sem tendenciosidade, com alta 
precisao e exatidao a baixo custo, foram feitos estudos no sentido de 
selecionar, entre os tamanhos dos segmentos estudados, o mais apropri_a 
do para ser utilizado como unidade amostra‘1 e ao mesmo tempo determinar 
o menor numero destes segmentos para compor a amostra. 

3.5.4 - SELECAO 00 TAMANHO DO SEGMENTO E PETERMINAQAQ DO Not™ DE UNI 
DADES POP AMOSTRA. 


0 sucesso do uso de esquema de amostragem para levan. 
tamento de areas agrTcolas nao depende somente da escolha do model o es^ 
tatTstico a ser utilizado. Dois fatores devem ser considerados num pU 
nejamento por amostragem: o tamanho do segmento (unidade amostra!) e o 
numero destas unidades que irao constituir a amostra. 

De urn modo geral , a escolha do tamanho do segmento e 
da amostra a ser utilizada em qualquer sistema de amostragem e baseada 
em duas premissas: estatTstica e economica. Do ponto de vista estate 
tico, o tamanho do segmento e da amostra devem ser suficientemente grari 
des para que as caracterTsticas das amostras representem as caracterTs 
ticas dos parametros populacionais correspondentes. Do ponto de vista 
economico, deve-se escolher urn tamanho de semento e de amostra, de tal 
modo que minimize o custo e o tempo dispendidos durante a sua execugao. 

Por outro lado, como esta pesquisa trata da obtengao 
de urn sistema de amostragem, onde uma das variaveis e a cl assificagao 
dos dados do MSS do LANDSAT, el a so se torna viavel quando a estimate 
va da area ocupada com trigo, obtida por este processo, for igual ou 
mais exata do que aquela obtida com o sistema 1-100. 
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Pel a comparagao da area ocupada com trigo, obtida atra. 
ves do slstema 1-100, utilizando-se dados do MSS do LANDSAT (30.447 ha), 
com a mesma area determinada pel a interpretagao visual de fotografias 
aeroas (28.314 Ha), verificou-se que essa diferenga em termos percen. 
tuais foi de + 7,5355. 

Baseado nestes conhecimentos, procurou-se nesta etapa 
fazer uni estudo com o objetivo de minimi zar o custo e o tempo dispendi^ 
do durante os levantamentos por amostragem, quando se utiliza a estima 
tiva pela regressao. Para o atingimento deste objetivo foram testados 
como unidades amostrais cinco tamanhos de segmentos 10, 20, 30, 40 e 60 
km 2 . Cada tamanho de segmento foi considerado como urn piano amostral . 
Procurou-se estudar cada piano amostral variando o tamanho da amostra 
(i.e., diferentes proporgoes de area amostrada) e para cada tamanho de 
amostra, 20 repetigoes foram extrandas da populagao, utilizando-se amos. 
tragem aleatoria simples sem repos igao. A Figura 3.2 mostra esquematica 
mente como foram obtidas as areas tritTcolas utilizadas na estimativa 
pela regressao para selecionar o tamanho da unidade amostral e o tamanho 
otimo da amostra. Porem, a utilizagao da estimativa pela regressao to£ 
na-se invalida quando e usada uma amostra com 2 unidades (n=2). Neste ca 
so, r 2 =l e, consequentemente, a variancia da estimativa assume o valor 
zero para qualquer que seja a area estimada. Desta forma, para que os re^ 
sultados nao fossem afetados por este fator, o menor numero de unidades 
tomadas para cada piano amostral foi 4. 

Uma vez obtidos os resultados das estimativas das areas 
ocupadas com trigo para as 20 repetigoes de cada tamanho de amostra, pr£ 
curou-se determinar, em cada piano amostral, a menor amostra que est£ 
mou a area desta cultura com boa precisao e exatidao. A determinagao de£ 
tas amostras foi feita com base em quatro criterios adotados: 
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Criterio 1 


Criterio 2 


Criterio 3 


Criterio 4 


- quc a diferenga relativa nao fosse maior que - 7 , 53% 
para a area obtida com a estimativa pel a regressao, 
diferenga esta nao maior do que aquela registrada 
pel a interpretagao digital; 

- que a media da populagao das diferengas entre as 
areas estimadas pelo sistema amostral e aquela obtida 
atraves de fotografias aereas fosse zero (a=Q,05); 

- que a estimativa obtida pelo sistema amostral nao fos_ 
se tendenciosa; 

- que o coeficiente de variagao nao fosse superior a 
5 %> ou seja, que o dcsvio padrao nao ultrapassasse 5 % 
da referida estimativa. 
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Fig. 3.2 - Esquema do processo de obtengao das areas tritTcolas. 






flnidos como: 


~ 41 *» 


A diferenga relativa e o coeficiente de variagao sao d£ 


1) - Diferenga Relativa (OR) 


DR (ha) «f R - 28, 314 ha 


ou em percentagem 


DR (55) = 


Y r - 28.314 ha 
~ 28. 314 ha 


X 100, 


onde 20.314 ha correspondente a area ocupada com trigo, obtida de foto 
grafias aereas. 


2) - Coeficiente de Variagao (CV) 


C.V. (55) 


Y R 


X 100. 


Os resultados das estimativas das areas ocupadas com 
trigo, obtidas de cada tamanho de amostra, foram plotados nun? conjunto 
de coordenadas cartesianas. Desta forma, para cada piano amostral tes_ 
tado obteve-se urn grafico, onde no eixo das abscissas variou-se o ta 
manho da amostra e no eixo das ordenadas, os valores das areas estima. 
das. Sobre cada grafico, tres retas paralelas ao eixo das abscissas f£ 
ram tragadas. Tais retas sao: a) uma reta central, a qual representou 
a area ocupada com trigo, obtida atraves de fotografias aereas e; b) 
duas retas que representaram as areas superestimadas e subestimadas com 
diferengas de - 7,5355. Atraves destes graficos e utilizando o criteri£ 
1, foram excluTdos todos os tamanhos de amostras que forneceram ; pelo 
menos uma estimativa, dentre as 20 repetigoes, com uma diferenga maior 
que 7,5355, quando comparado com o parametro populacional , Baseado no 



crl terlo-1 , fez-se a selegao do manor tamonho da amostra do cada piano 
amostral, Posteriormente, para verificar so a media das diferengas re, 
lativas, obtidas de 20 repetigoes para cada tamanho do amostra selccio 
nada, foi Iguala zero, conforme o criterio-2, nplicou-se o teste t 
(a c 0,05). Satisfeito o segundo criterio, foram tomadas as estimativas 
das amostras selocionadas e submetidas a uma anallse de tendenciosida, 
de de acordo com o criterio-3. A analise de tendenciosidade foi verifi 
cada atraves do teste do sinal com a » 0,05, (t-loel, 1966). Finalmente, 
determinou-se para cada estimativa o coeficiente de variagao, a fim de 
verificar, com base no criterio-4, se os valores destes coeficientes 
nao ultrapassavam a 5% . Admitiu-se como o menor tamanho da amostra aque 
le que satisfez os quatro criterios sinuiltaneamente. 

Para quantificar a eficiencia da estimativa em que fo 
ram utilizados os dados auxiliares do satelite LANDSAT,procurou~se com 
pararasareas tritTcolas, obtidas atraves da estimativa pel a regressao, 
com as areas estimadas usando o metodo de expansao direta (somente da_ 
dos de aeronave). Em ambos os metodos foram utilizadas as mesmas unidia 
des amostrais para a aquisigao dos dados necessarios a estimativa da 
area. A expressao matematica do metodo de expansao direta usada neste 
estudo e: 

n 

^ED = y-j/ n ) 9 


onde: 


- area ocupada com trigo, obtida dos dados de aeronave para a 
unidade amostrada i, i = 1 a n. 


N - numero de unidades amostrais da populagao. 
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formul a : 


A variancia da area estimada foi obtida pela seguinte 



N - n 
N 


{ (yj - y) 2 - 


A Eficiencia Relativa (ER) da estimativa pela regre£ 
sao comparada a expansao direta foi medida atraves da razao das respe£ 
tivas variancias, ou seja: 

E.R, B V(V ED ) / V (? R ). 

A necessidade de repetigao de amostra e a precisao nos 
calculos matenmticos foram considerados fatores importantes na conipa 
ragao dos resultados obtidos nas etapas de selegao do tamanho do seg 
mento e na determinagao do numero de unidades por amostra, Como a es 
timativa da area de trigo atraves do metodo de estimativa pela regre£ 
sao exige calculos relativamente complexes, todos os calculos matema 
ticos, bem como a extragao de amostra, foram feitos pelo copiputador, 
confornie o programa apresentado no Apendice B. 



CAPITULO 4 


RESULTADOS E DISCUSSftO 


4.1 - RESULTADOS ALCANQADOS COM DADOS DE AERONAVE E DO LANDSAT 


Os temas intorpretados a parti r das fotografias aereas 
foram: trigo, solo preparado, pastagem, mata e o tema outros (Tabela 
4.1), 


TABELA 4,1 

TEMAS INTERPRETADOS A PARTI R DE FOTOGRAFIAS AEREAS, DIVIDIDOS 

EM CLASSES E SUBCLASSES 




1 

TEMA 

STMBOLO 

CLASSE 

SUB CLASSE 



T1 



T2 

TRIGO 

T 

T3 



T4 



PI 

PASTAGEM 

P 

P2 

SOLO PREPARADO 

SP 

- 

MATA 

M 

- 

OUTROS 

0 

- 
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0 tenia trigo apresentou diferentes padroes de tonal idades 
que variaram de magenta escura a clara. Devido a essa diferenga no padrao 
tonal, o tema trigo foi dividido ein 4 subclasses denominadas Tl, T2, T3 
e T4. Posteriormente, com a observagao de campo,verificou-se que a varia 
gao no padrao de tonal idade foi, principalmente, devido a diferentes varia 
dades de trigo e a diferentes estaaios de desenvol vimento. 

0 tema pastagem foi considerado de forma general izada iri 
cluindo pastagem natural, cultivada e areas eni pousio. A niaior parte da 
pastagem apresentou-se bastante degradada na epoca do aerolevantamento. 
Nos locais onde houve condi goes mais favoraveis de umidade, constatou-se 
a presenga de pastagem mais exuberante, Tal fato acarrotou o aparecimento 
de dois padroes de tonal idades quando observados atraves de fotografias 
aereas: 1) magenta bem Clara mosqueada de verde, que e caracterTstica de 
pastagem degradada ou areas em pousio, e, 2) magenta cl ana mosqueada por 
outra magenta mais escura, semelhante a urn dos padroes tonais de trigo 
para as areas onde a pastagem apresentou maior vigor vegetativo. Por e_s 
ta razao o tema pastagem foi dividido em 2 subclasses denominadas PI e 
P2. 


No caso dos tal hoes de pastagem que apresentaram urn pa_ 
drao de tona'lidade semelhante ao do trigo, a separagao visual dos dois te 
mas (trigo e pastagem) foi possTvel utilizando as caracteristicas espaci_ 
ais, texturais e informagoes complementares como: curvas de nTvel, car 
readores, presenga de gado e localizagao da area problema com relagao aos 
a! vos vizinhos 


0 tema outros incluiu cidade.estradas e rios presentes 
na area de estudo. Este tema junto com os temas mata e solo preparado, por 
apresentar padroes de tonal idade distintos do tema trigo, nao foram di vj_ 
didos em subclasses. 

0 tema mata teve como fator preponderante para asuaideji 
tificagao a textura rugosa caracterTstica. A principal razao para sua 
interpretagao foi servir como urn dos pontos de referencia para a local^ 
zagao da area de estudo nos dados do LANDSAT. 
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Durante a fase de interpretagao digital, em que foi uti_ 
lizado o metodo HTbrido M-2, as diferentes feigoes da area de estudo 
foram agrupadas, pela abordagem K-MEDIAS ("clustering"), em oito cla£ 
ses espectrais homogeneas. Essas classses foram assim distribuTdas:qua 
tro classes espectrais para o tema trigo, duas classes para o tenia pas_ 
tagem, uma classe espectral para o tema mata e uma classe para o solo 
preparado. 


A separagao dos temas trigo e pastagem em quatro e duas 
classes espectrais, respectivamente, pela abordagem K-MEDIAS, serviu pa_ 
ra mostrar que as varias tonal idades do trigo e pastagem manifestadas 
nas fotos aereas tambem foram constatadas -pelo sensor MSS do LANDSAT d£ 
vido a mudanga de nTvel de cinza. Isto levou a conclusao de que as clas. 
ses trigo e pastagem nao apresentaram urn uni co comportamento espectral 
de forma gaussiana, o que e pre-requisito para o uso do classificador 
MAXVER e, portanto, foram di vi di das em subclasses espectrais. 

0 tema "outros", constituTdo em sua maior parte pela 
cidade de Cruz Alta, foi agrupado nos temas pastagem e solo preparado 
por apresentar muito solo exposto e algumas areas com vegetagao;porta_n 
to, espectralmente semelhante a esses temas. 

Embora as cuK'uras de cevada e centeio influenciem os 
resultados da classificagao digital da cultura do trigo, no municTpio 
de Cruz Alta, segundo a Fundagao Instituto Brasileiro de Geografia e Es_ 
tatTstica - IBGE (1979) no ano de 1979 nao houve plantio de centeio e 
a area plantada com cevada foi apenas 104 ha, o que representa, em ter 
mos percentuais, 0,03% da area daquele municTpio. 

Nas Tabelas 4.2 a 4.6 sao mostrados, para cada segmen_ 
to estudado (10, 20, 30, 40 e 60 km 2 ), os resultados comparatives das 
areas ocupadas com trigo, obtidas a partir da interpretagao das fotos 
aereas e da classificagao digital utilizando dados do MSS do LANDSAT. 
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TABELA 4,2 

AREAS OCUPADAS COM A CULTURA DO TRIGO, OBTIDAS ATRAVES DE DADOS DO 


I.ANDSAT E DE FOTOS AEREAS, PARA OS SEGMENTOS DE 10 kin 2 


IBSBSBSI 

LANDSAT 

FOTO AEREA 

<160,-14 

435,65 

4 G2 ,63 

460,14 

606,76 

570,62 

525,83 

486,4?. 

604,35 

568,53 

395,05 

396,64 

716,04 

650,35 

634,23 

644,08 

574 ,80 

551,51 

478,66 

487,82 

469,70 

402,81 

286,36 

325,63 

639,30 

498,66 

368,47 

340,70 

467,31 

467,01 

208,72 

200,06 

444,32 

553,60 

416,25 

452,98 

565,25 

542,26 

473,26 

426,70 

272,32 

235.89 

153,18 

120,34 

213,80 

239,18 

294,12 

338,01 

587,64 

570,62 

553,31 

550,92 

443,42 

441,33 

392,96 

365,19 

93,46 

65,69 

482,54 

487,32 

282,18 

268,44 

48,67 

34,94 

260,68 

278,59 

212,90 

194,39 

44,19 ■ 

20,60 

303,08 

272,97 
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TABELA 4.3 

*EAS OCUPADAS COM A CULTURA DO TRIGO. OBTIDAS ATRAVES DE DADOS DO 


LANDSAT E DE FOTOS AEREAS. PARA OS SEGMENTOS DE 2Q km 2 


N9 DO 
SEGMENTO 





AREA DE 

TRIGO (ha) 

N9 DO 

LANDSAT 

FOTO AEREA 

SEGMENTO 

1120,05 

935,22 

19 

1076,45 

950,74 

20 

1076,45 

970,45 

21 

923,27 

893,41 

22 

1132,59 

1059,04 

23 

999,41 

965,37 

24 

1314,73 

1132,59 

25 

1114,08 

885,65 

26 

976,42 

909,24 

27 

1350,26 

1294,43 

' 28 

1053,46 

1038,83 

29 

756,05 

738,44 

30 

1077,95 

990,16 

31 

1298,91 

1191,12 

32 

1108,10 

1014,04 

33 

1007,78 

839,36 

34 

676,06 

667,07 

35 

860,56 

1006,58 

36 



476,86 

653.63 

850.41 
1043,60 

427,60 

507,91 

635.42 
866,53 

590.63 
1140,95 

836,38 

576,00 

977,02 

570,62 

219,17 

330,85 

473,57 

347,27 


286,06 

546,74 

707.68 
968,96 
356,26 
577,19 

644,42 
786,81 
547,93 

1121,54 

806,51 

533,00 

962.69 
624,07 
269,33 
303,38 
472,98 
294,12 
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TAB EL A 4.4 

AREAS OCUPADAS COM A CULTURA DO TRIGO, OBTIOAS ATRAVES DE DADOS DO 


LANDSAT E DE FOTOS AEREAS, PARA OS SEGMENTOS DE 30 kin 2 


N9 DO 
SEGMENTO 


AREA DE 

TRIGO (ha) 

LANDSAT 

FOTO AEREA 

1699,33 

mmm 

1573,62 

WBm 

1530,02 


1525,25 

WfflM 

1834,90 

1552,13 

1570,34 

1375,35 

.1925,07 

1845,94 

1234,71 

1225,75 

1789,21 

1667,68 

1695,75 

1527,64 

1475,08 

1306,38 

1069,64 

1206,64 


N9 DO 
SEGMENTO 


AREA DE 

TRIGO (ha) 

LANDSAT 

FOTO AEREA 

777,55 

566,74 

1203,36 

973,73 

1318,02 

1204,85 

661,10 

697,53 

986,28 

944,77 

1106,31 

1034,05 

1584,37 

1562,87 

968,96 

918,20 

1352,96 

1388,79 

413,86 

467,31 

591 ,53 

' 581,97 

560,17 

488,51 
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TABELA 4.5 


AREAS OCUPADAS COM A CULTIJRA DO TRIGO, OBTIDAS ATRAVES DE DADOS DO 


LANDSAT E DE FOTOS AEREAS, PARA OS SEGMENTOS DE 40 km 2 


N9 DO 
SEGMENTO 



AREA E 

IE TRIGO (ha) 

N9 DO 

Area de 

TRIGO (ha 

LANDSAT 

FOTO AEREA 

SEGMENTO -| 

LANDSAT 

FOTO AEREA 

2434,78 

2067,80 

10 

2273,54 

2190,22 

1554,81 

1276,22 

11 

2051 ,38 

1808,32 

1612,44 

1610,77 

' . ■ 

1471,80 

1424,92 

2190,52 

1836,39 


2186,05 

2095,87 

1952,53 

1737,85 

14 

1103,62 

1023,30 

1437,16 

1410,88 

15 

1309,95 

1279,50 

2052,87 

1879,69 

16 

1755,47 

1703,81 

1958,50 

1721 ,72 

17 

1368,49 

1583,67 

809,80 

817,27 

18 

923,27 

847,13 


TABELA 4.6 

AREAS OCUPADAS COM A CULTURA DO TRIGO, OBTIDAS ATRAVES DE DADOS DO 


LANDSAT E DE FOTOS AEREAS, PARA OS SEGMENTOS DE 60 km 2 


AREA I 

IE TRIGO (ha) 

N9 DO 

AREA 

DE TRIGO (ha! 

’LANDSAT 

FOTO AEREA 

SEGMENTO ’ 

LANDSAT 

FOTO AEREA 

3534,23 

2928,10 

7 

2793,10 

2511,23 

3143,96 

2798,78 

8 

1730,38 

1904,17 

3455,10 

3312,37 

9 

2339,23 

2333,56 

2759,95 

2677,55 

10 

1520,17 

1501,36 

'2566,76 

2234,42 

11 

2175,90 

2144,84 

2899,11 

2501,37 

12 

1529,13 

1406,70 
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Muito embora o metodo HTbrido M-2 tenha side o melhor 
entre os metodos estudados (Lima et alii, 1982), houve uma certa con 
fusao da classo trigo com a classe pastagem e, em menor intensidade, 
com a classe mata (Tabela 4.7). Isto se deve ao fato de que, na ep£ 
ca da tomada dos dados do LANDSAT, alguns talhoes de pastagem mais 
exuberantes tiveram resposta espectral semelhante a do trigo. Para o 
caso da classe mata, notou-se que, na maioria das vezes onde havia 
confusao com o trigo, a vegetagao era inuito esparsa, , com predominan 
cia de urn estrato inferior graminoide, o que provavelmente causou uma 
semelhanga de resposta espectrais. 

TABELA 4.7 

M ATRIZ DA CLASSIFICAQflO 


IBM 

as 

1 

2 

3 

n 

5 

6 

a 

8 

1 . TRIGO 1 

ESI 

100,0 

0,0 

0,0 

0,0 

o 

*x 

o 

0,0 

0,0 

0,0 

2. TRIGO 2 

0,0 

0,0 

99,0 

0,8 

0,0 

0,0 

0,0 

0,0 

0,0 

3. TRIGO 3 

0,0 

3,1 

1,6 

93,7 

0,0 

0,0 

0,0 

1,6 

0,0 

4. MATA 

0,0 

0,0 

0,0 

0,0 

96,1 

0,0 

0,0 

2,6 

1,3 

5. SOLON U 

0,0 

0,0 

EO 

0,0 

0,0 

97,9 

0,0 

eh 

0,0 

6; PASTO 1 

1,0 

0,0 

0,0 

0,0 

0,0 

0,0 

86,5 

12,5 

0,0 

7. PASTO 2 

0,0 

0,4 

0,4 

1,1 

0,6 

0,0 

8,5 

88,6 

0,4 

8. TRIGO 4 

0,0 

0,0 

0,0 

0,0 

0,0 

0,0 

0,0 

6,2 

93,7 


* - N = pontos da area amostrada que nao foram claeaificadas em nenhuin 
doe temasj 

- os numeroe de 1 a 8 na horizontal da matriz correepondem aos te_ 
mas contidos na linha vertical ; 

- os numeroe aontidos na diagonal principal indicam a peroentagem 
da ciaseificagdo correta de cada terna contido na matriz ; 

- oe demais numeroe da matriz indicam o grau de confusao entre os 
temas. 
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Conforme se observa nas Tabelas 4.2 a 4.6, de um modo gie 
ral as areas estimadas com a cultura do trigo, atraves do 1-100, util 1za_n 
do os dados do MSS do LANDSAT foram superestimadas quando comparadas com 
as mesmas areas obtidas atraves da interpretagao visual das fotografias 
aereas. A causa desta superestimativa foi atribuTda ao problema de reso_ 
lugao do si sterna sensor do LANDSAT, bem como a confusao de respostas es_ 
pectrais dos temas mata e pastagem com o tenia trigo. Em alguns segmentos 
as areas tritTcolas. foram subestimadas. Isto pode ser explicado pelo fa. 
to de que a pastagem nestes segmentos encontrava-se em estado bastante 
degradado, nao sendo espectralmente confundida com o trigo (baixo erro 
de inclusao). Por outro lado, a baixa densidade da cultura do trigo em 
alguns talhoes destes segmentos, ocasionada por falha no plantio, forn£ 
ceram uma resposta espectral semelhante a resposta espectral de solo pre. 
parado, causando assim o erro de omissao durante a classificagao digital. 
Os altos erros de omissao nao foram compensados pelos erros de inclusao, 
consequentemente, houve uma subestimativa da area de trigo nestes segmeri 
tos. 


4.2 - ESTIMATIVA DA AREA OCUPADA COM TRIGO PELO SISTEMA DE AMOSTRAGEM 
4.2.1 - ANflLlSE DA CORRELAQAO 


A Tabela 4.8 mostra os coeficientes de corregao para os 
segmentos estudados. 
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TABELA 4,8 

COEFICI ENTES DE CORRELATE) ENTRE AREAS DE TRIGO ESTIMADAS USANDO DADOS 

DO LANDSAT E FOTOGRAFIAS AEREAS 


TAMANHO DO 
SEGMENTO 
(km 2 ) 

N9 DE OBSERVA 
gOES 

COEFICIENTE DE CORRELAgAO 

SENTIDO 

HORIZONTAL 

VERTICAL 

DIAGONAL 

60 

12 

0,97 

0,98 


40 

18 

0,97 

0,98 

0,98* 

30 

24 

0,97 

0,98 

0,96 

20 

36 

0,97 

0,98 

0,96 

10 

72 

0,96 

0,96 

0,96 


* E8tc coeficiente foi calculado ubi'lizandjo eegmentoa da 2 x 2 wiidadaa 
baaica8 da 10 ton 2 . 

De acordo com os resultados observados na Tabela 4.8,n£ 
ta-se que, independente do tamanho dos segnientos e do sentido, todos os 
coeficientes de correlagao foram altamente sigm'ficativos ({a = 0,05) e 
tiveram valores iguais ou acima de 0,96. Isto niostra que, para os tama_ 
nhos de segmentos estudados, ha uma alta correlagao entre a area estima 
da com dados do LANDSAT e a area obtida atraves das fotografias aereas. 

4.2.2 - ANALISE DE REGRESSAO SIMPLES 


Devido a alta correlagao (r .>0,96) obtida, a analise de 
regressao simples foi feita apenas para os dados dos segnientos de 10, 
20 e 30 km 2 tornados no sentido horizontal e para os segnientos de 40 e 
60 km 2 foram utilizados os resultados obtidos somando 2 x 2 e 3 x3 unf 
dades de 10 km 2 , respectivamente. Os graficos das Figuras 4.1 a 4.5 mos 
tram a correlagao existente entre as areas estimadas com a cultura do 
trigo, obtidas atraves de dados do MSS do LANDSAT e dados de fotografias 
aereas. 
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Lima vez determinadas as equagoes das retas que ropresen^ 
taram a relagao entre as variaveis, testaram-se os coeficientes de inter 
sc?ao das retas para veriflcar se eram diferentes de zero, Para isto, 
aplicou-se o teste t (a = 0,05), cujo resultado mostrou que os coefi^ 
cientes de Intersegao das retas para os cinco pianos aniostrais diferj[ 
ram estatistlcamente de zero. 



AREA ESTIMADA COM DADOS DE FOTO AEREA (x TOO ha) 
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ORIGINAL PAGE 19 
OF POOR QUALITY 



Fig. 4.1 - Regressao entre as estimativas das areas de trigo, obtidas 
atraves de dados do LANDSAT e dados de aeronave, para o ta 
manho de segmento de 10 km 2 . ~ 
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original page rsr 

OF POOR QUALITY 



Fig. 4.2 - Regressao entre as estimativas das areas de trigo, obtidas 
atraves de dados do LANDSAT e dados de aeronave, para o ta 
manho de segmento de 20 km?. 




AREA 0 
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ORIGINAL PAGE 10 
OF POOR QUALITY 



Fig. 4.3 - Regressao entre as estimativas das areas de trigo, obtidas 
atraves de dados do LANDSAT e dados de aeronave, para o ta 
manho de segmento de 30 km 2 . 



S1EA ESTIMADA COM DADOS DE FOTO AEREA {x ICO ha} 
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original PW3EI3 



Fig. 4.4 - Regressao entre as estimativas das areas de trigo, obtidas 
atraves de dados do LANDSAT e com dados de aeronave, para 
o tamanho de segmento de 40 km 2 . 






SREA ESTIMADA COM DADOS DE FQTO AEREA (x 100 ha) 
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ORIGINAL PAGE {S 
OF POOR QUALITY 



Area estimada com dados do landsat (x too ha) 


Fi,0. 4,5 - Regressao entre as estimativas das areas de trigo, obtidas 
atraves de dados do LANDSAT e com dados de aoronave, para 
o tamapho de segmento de 60 km 2 . 
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Verificando que os resultados obtidos do LAMDSAT e das 
fotos aereas apresentaram uma alta correlagao (r r > 0,96) e as equa 
goes lineares nao passaram pel a origem dos eixos de acordo com indi_ 
'cagoes de Cochram (1965), optou-se pelo metodo de regressao para a 
estimativa da area de trlgo na regiao de estudo. 

4.2.3 - SELECWO DO TAMANHO DO SEGMENT'D E DETERMINAgAO DO NOMERO DE 
UNIDADES PQR AMOSTRA 

As Figuras 4.6 a 4.10 contem as estimativas das areas 
ocupadas com trigo, obtidas pelo metodo de regressao para os cinco 
pianos amostrais e variando-se o tamanho da amostra. Para eada tatirn 
nho de amostra foram obtidas estimativas de 20 repetigoes. 

Com base no primeiro criterio adotado na CapTtulo 3, 
excluiram-se todos os tamanhos de amostras que forneceram pelo me nos 
uma estimativa com diferenga superior a 7,53$. Nas Figuras 4.6 a 4.10 
pode-se ver que o me nor tamanho de amostra (para cada piano amostra!)' 
que satisfez este criterio foi composto de 10, 8, 8, 6 e 5 unidades 
para os tamanhos de segmentos de 10, 20, 30, 40 e 60 km 2 . Em termos 
percentuais de area amostrada isto equivale a 13,90%; 22,22%; 33,33%; 
33,33% e 41,67% respectivamente. 

Mas Tabelas 4,9 a 4.13 estao os resultados das areas 
triticolas, obtidas atraves da estimativa pela regressao, para os fra 
manhos de amostras selecionadas. Mostra-se, tainbem, nestas tabelas o 
valor da diferenga relativa e do coeficiente de variagao para a area 
estimada em cada repetigao. 



est:*at:va ca a'p.ea TaiTfcoi_A cr. io* 
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ORIGINAL PAG£ *3 
OF POOR QUALITY 



Fig. 4.6 - Estimativa da area triticola, obtida atraves da estimativa 
pela regressao, usando segmentos de 10 km 2 e variando o 
tamanho da amostra. 




ESTIMATIVA DA AREA TRITlCOLA 
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ORIGINAL PAGE 13 
OF POOR QUALITY 



Fig. 4.7 ~ Estimativa da area tritTcola, obtida atraves da estimativa 
pel a regressao, usando segmentos de 20 km 2 e variando o 
tamanho da amostra. 





ESTJf.'ATIVA DA AREA TR IT ICO LA 
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ORIGINAL PAGE IS 
OF POOR QUALITY 



Fig. 4.8 - Estimate va da area tritTcola, obtida atraves da estimativa 
pel a regressao, usando segmentos de 30 km 2 e variando o 
tamanho da amostra. 






ORIGINAL PAGE J0 
OF POOR QUALITY 
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Fig. 4.9 - Estimativa da area tritTcola, obtida atraves da estimativa 
pel a regressao, usando segmentos de 40 km 2 e variando o 
tamanho da amostra. 
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ORIGINAL PAGE 19 
OF POOR QUALITY 



Fig. 4.10 - Estimativa da area tritTcola, obtida atraves da estimativa 
pel a regressao, usando segmentos de 60 km 2 e variando o 
tamanho da amostra. 
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TABULA 4.9 


RESULTADOS OBTIDOS ATRAVE.S DA ESTIMATIVA PELA REGRESSAO, PARA 0 SEGHENTO 
DE TO km ? -, USANDO AMOSTRAS COM 10 UNIDADES 


N9 DE 

REPETICAO 

AREA ESTIMADA 
? R (ha) 

DIFERENfA RELATIVA 

VARIANCIA 

vtV 

c.v. 

(%) 

(ha). 

C«1 

1 

28.706,75 

392,75 

1 ,39 

943.101,31 

3,4 

2 

28.963,22 

649,22 

2,29 

360.948,87 

2,1 

3 

27.808,36 

-505,64 

-1,78 

311.541,23 

2,0 

4 

29,259,53 

945,53 

3,34 

437.509,15 

2,3 

5 

28.316,01 

2,01 

0,01 

1 .013.768,71 

3,5 

6 

27.384,65 

-929,35 

-3,28 

800.655,31 

3,3 

7 

29.082,70 

768,70 

2,71 

566.055,18 

2,6 

8 

29.483,67 

1169,67 

4,13 

602.561,09 

2,6 

9 

28,576,66 

262,66 

0,93 

433.545,47 

2,3 

10 

28.447,07 

133,07 

0,47 

439.066,36 

2,3 

11 

28.180,59 

-133,41 

-0,47 

1 .010,956,75 

3,6 

12 

26.228,49 

-2085,51 

-7,36 

1.174.274,80 

msm 

13 

29.542,21 

1228,21 

4,34 

1.789.146,87 


14 

29.287,94 

973,94 

3,44 

1.560.397,32 


15 

28.716,81 

402,81 

1,42 

246.227,12 

H 

16 

28,129,50. 

184,50 

-Q,65 

1 .163.725,68 

3,8 

17 

27.609,15 

-704,85 

-2,49 

724.708,85 

3,1 

18 

29.610,69 

1296,69 

4,58 

649.077,00 

2,7 

19 

28.272,21 

-41,79 

-0,15 

1.024.127,30 

3,6 

20 

28.515,89 

201,89 

0,71 

1.044.772,71 

3,6 

MEDIA (X) 
VARIANCIA 

r \l fV\ 

28.506,10 
691 .475,40 

650,61 

2,30 

- 

3,1 


TESTE T 
R.M.Q.E. (ha) 


2,92 


1 ,03n .s . 
832,95 


2,29 
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TABELA 4.10 


RESULTADOS OBTIDOS ATRAVES DA ESTIMATIVA PELA REGRESSAO, PARA 0 SEGMENTO 


DE 20 km z , USANDO AMOSTRAS COM 8 UMIDADES 



C.V. (%\ 
TESTE T 


AREA ESTIMADA 

DIFERENQA RELATIVA 

?„ (ha) 

(ha) 

(*> 

28.875,41 

561,41 

1,98 

29.826,20 

1512,20 

5,34 

27.587,45 

-726,55 

-2,57 

29.151 ,98 

837,98 

2,96 

28.299,64 

-14,36 

-0,05 

26.592,84 

-1721,16 

-6,07 

29.453,62 

1139,62 

4,02 

27.958,59 

-355,41 

-1,25 

28.754,32 

440,32 

1,56 

27.321 ,28 

-992,72 

-3,51 

29.247,14 

933,14 

3,29 

27.661,01 

-652,99 

-2,31 

26.828,52 

-1485,48 

-5,25 

28.214,15 

-99,85 

-0,35 

27.010,99 

-1303,01 

-4,60 

28.277,45 

-36,55 

-0,13 

28.101,66 

-212,34 

-0,75 

28.590,04 

276,04 

0,97 

30.196,82 

1882,82 

6,65 

28.086,46 

-227,54 

-0,80 

} 28.301 ,78 

770,37 

2,70 

A 964.579,34 

- 

- 

3,47 

- 

- 

- 

0,06n.s . 

- 

{ha I- 

957,34 

3,40 




398.747,86 
1 .077.848,66 
995.891,81 
1 .178.863,43 
1 .262.210,63 
676.031,77 

1.054.323.05 

670.903.69 
1 .885.307,96 

586.200,67 
1 .447.373,42 
368.902,99 

570.928.69 
532.665,56 

1.047.512.06 
478.294,11 
664.734,14 
459.600,24 
868.416,72 
484.065,75 












TABELA 4.11 

RESULTADOS OBTIDOS ATRAVES DA ESTIMATIVA PELA REGRESSAO, PARA 0 SE6MENT0 
DE 30 km?, USANDO AMOSTRAS COM 8 UNIDADES 


DO AREA ESTIMADA 


REPETIQAO 


1 




MtDIA 00. 
VARIANCIA 2 

c.v. t%\ 

TESTE T 
R.M.Q.E. (ha) 


? R (ha). 


28.080,35 

27.268.18 

27.945.88 

29.239.42 

29.244.63 
28.971,12 
29.309,83 

28.646.64 
28.450,67 
28.411,08 

28.594.32 

28.209.88 
28.031,59 

28.283.32 
28.357,62 
28.526,30 

28.296.42 
29.064,15 

28.744.18 
28.378,56 


28.502,71 

241.248,65 

1,72 


DIFERENQA RELATIVA 

VARIANCIA 

C.V. 

(ha) 

(*) 

Wr> 

(*) 

-233,65 

-0,83 

238.538,20 

1,7 

-1045,82 

-3,70 

578.416,06 

2,8 

-368,12 

-1,30 

611.190,37 

2,8 

925,42 

3,30 

59 .230,74 

0,8 

930,63 

3,30 

741 .227,49 

2,9 

657,12 

2,30 

284.251,20 

1,8 

995,83 

332,64 

3,50 

1,20 

575.117,93 

170.687,01 

2,6 

1,4 

136,67 

0,50 

386.614,85 

2,2 

97,08 

0,30 

331 .083,47 

2,0 

280,32 

0,90 

111.214,70 

1,2 

-104,12 

-0,40 

193.353,12 

1,6 

-282,41 

-0,90 

287.326,96 

1,9 

-30,68 

-0,10 

467.394,06 

2,4 

43,62 

0,20 

704.255,10 

2,9 

212,30 

0,70 

982.042,05 

3,5 

-17,58 

-0,10 

882.770,34 

3,3 

750,15 

2,60 

812.477,13 

3,1 

430,18 

1,50 

470.198,92 

2,4 

64,56 

0,20 

283.225,37 

1,9 

396,94 

1,40 

*• 

2,3 

- 

- 

- 

- 


- 

- 

** 

l,67n.s. 

- 

- 

- 

526,20 

1,80 

- 

- 
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TABELA 4.12 


RESULTADOS OBTIDOS ATRAVES DA ESTIMATIVA PELA REGRESSAO.PARA 0 SEGMENTO 


DE 40 km 2 , >USANDO AMOSTRAS COM 6 UNIDADES 


N9 DE 
REPETig/tt) 



'Area estimada 

DIFERENCA relativa 

VARIANCIA 

? R Cha) 

(ha) 

(*) 

v<r R ) 

29.076,46 

762,46 


480.220,12 

27.969,64 

-344,36 


1 .371 .886,69 

27.286,79 

-1027,21 


365.045,15 

27.327,88 

-986,12 

-3,50 

705.299,44 

29.593,39 

1279,39 

4,50 

746.008,96 

30.183,27 

1869,27 

6,60 

627.255,41 

29.217,57 

903,57 

3,20 

218.033,03 

28.781,29 

467,29 

1,60 

303.271,27 

28.925,98 

611,98 

2,20 

1 .041 .660,17 

28.763,43 

449,43 

1,60 

1 .033.297,65 

28.549,43 

235,43 

0,80 

353.526,31 

28.524,07 

210,07 

0,70 

649.903,19 

27.035,87 

-1278,13 

-4,50 

337.768,43 

27.745,68 

-568,32 

o 

o 

CM 

K 

384.183,83 

28,003,86 

-310,14 

-1,10 

225.370,28 

29.144,68 

830,68 

2,90 

746.556,46 

28.351,36 

37,36 

0,10 

888.905,83 

27.230,81 

-108,19 

o 

CO 

CO 

1 

801 .099,45 

27.749,27 

-564,73 

-2,00 

425.839,80 

27.972,91 

-341,09 

-1,20 

483.812,40 


MEDIA C U 28.371,68 
VARIANCIA 692.781,15 
C.V. C %\ 2,90 
TESTE T 

R.M.Q.E. (ha) - 


708,01 


0,3Qn.s . 
834,33 
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TABELA 4.13 


RESULTAPOS OBTIDOS ATRAVES DA ESTIMATIVA PELA REGRESSAO.PARA 0 SEGMENTO 


DE 60 km 2 , USANDO AMOSTRAS COMO 5 UNIDADES 


de area estimada 



REPETigAO 


1 


27.491,47 

28.059.56 
27.951 ,11 

28.476.56 

28.884.31 

27.844.32 
28.185,72 
28.805,05 
28.130,46 
27.940,77 
29.513,20 
29.486,98 
29.049,66 
29.502,96 
28.101,75 
28.852,36 
28.968,82 
29.513,20 
28.962,55 
28.687,52 


MEDIA 28.620,41 

VARIANCIA 396.246,42 
C.V. {%} 2,20 

TESTE T 
R.M.Q.E. (ha) 


DIFERENCA RELATIVA 

(ha) 

w 

-822,53 

-2,90 

-254,44 

-0,90 

-362,89 

-1,28 

162,56 

0,57 

570,31 

2,01 

-469,68 

-1,66 

-128,28 

-0,45 

491,05 

1,73 

-183,54 

-0,65 

-373,23 

-1,32 

1199,20 

4,23 

1172,98 

4,14 

735,66 

2,60 

1188,96 

4,20 

-212,25 

-0,75 

538,36 

1,90 

654,82 

2,31 

1199,20 

4,23 

648,55 

2,29 

373,52 

1,32 


587,10 2,07 


2,18n.s. 

685,89 2,42 


VARIANCIA 

V(? R ) 


416.944.94 

182.971.98 
148.870,84 
428.837,75 

386.422.98 
476.131,21 
265.105,02 

638.259.99 

555.612.95 
494.485,12 
449.490,08 

136.222.99 
364.680,38 
560.852,44 
226.563,62 

36.545,64 

223.694,97 

449.490,08 

767.202,01 

272.647,07 
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Para a analise do segundo criterio (l.e., verificar sg 
a media da populagao das diferenqas relativas nao diferiu estatTstica 
mente de zero), aplicou-se o teste t, cujo resultado nao revelou dif <2 
renga significativa ao nTvel de 5$ para as estimativas provenientes das 
amostras selecionadas dos segmentos de 10, 20, 30 e 40 km 2 . Entre 
tanto, as estimativas obtidas pel as amostras com 5 unidades usandooseg 
mento de 60 km 2 nao satisfizeram o criterio proposto. Por esta razao, 
excluiu-se este tamanho de amostra, niesmo atendendo ao primeiro crite 
rio, e para analises posteriores, foram utilizados os resultados obti. 
dos com a amostra de 6 unidades que satisfizeram aos criterios 1 e 2 
(Tabela 4.14). Assim, a percentagem da area amostrada, usando este se£ 
mento, passou de 41 ,67ft para 50$, e os menores tamanhos de amostras que 
satisfizeram os dois primeiros criterios foram compostos de 10, 8,8,6 , 
e 6 unidades para os segmentos de 10, 20, 30, 40 e 60 km 2 , respective 
mente. 


Pnsteriormente, procurou-se verificar se esses tamanhos 
de amostras selecionadas forneceram estimativas tendenciosas (tercej. 
ro criterio). Sobre este aspecto, procurou-se comparar, atraves do tes. 
tes do sinal (a c 0,05), as estimativas da area tritTcola obtidas atra^ 
ves de cada amostra com o parametro populacional (28.314 ha). Os re. 
sultados do teste mostraram que nenhum dos tamanhos de amostras sele^ 
cionados forneceram estimativas tendenciosas. 

Por outro lado, nota-se pel as Tabelas 4.9 a 4.14 que as 
20 estimativas das areas ocupadas com trigo, obtidas atraves da estima 
tiva pel a regressao para cada amostra selecionada, foram bastante co]i 
sistentes, o que foi demonstrado pelos baixos valores dos coeficientes 
de variagao, entre 1,72 e 3,47$. 

Verificando ainda, pelas Tabelas 4.9 a 4.14, as colunas 
que consti tuem os valores dos Coeficientes de Variagao (CV), nota-se 
que nenhum tamanho de amostra selecionado forneceu uma estimativa de 
area com CV maior que 5$, o que satisfez plenamente o quarto criterio 
adotado. 


n 


TABELA 4,14 

RESULTADOS OBTIDOS ATRAVES DA ESTIMATIVA PELA REGRESSAO,PARA 0 SEGMENTO 


DE 60 km 2 , USANDO AHOSTRAS COM 6 UNIDADES 





MEDIA 

variancia 

C.V. l%l 
TESTE T 


AREA ESTIMADA 

DIFERENCA RELATIVA 

VARIANCIA 

? R Cha) 

(ha) 

(*) 


20.174,66 

-139,34 

-0,49 

395.912,74 

28.106,89 

-207,11 

-0,73 

316.724,47 

28.099,82 

-214,18 

-0,76 

447.380,92 

28.145,60 

-168,40 

-0,59 

493.006,01 

29.176,99 

862,99 

3,05 

433.906,74 

28.189,50 

-124,50 

-0,44 

392.504,32 

28.339,71 

25,71 

0,09 

222.778,53 

29.791,36 

1477,36 

5,20 

186.740,28 

27.423,56 

-890,44 

-3,14 

236.185,81 

28.795,41 

481,41 

1,70 

135.551 ,82 

27.979,57 

-334,43 

-1,18 

170.816,86 

28.348,21 

34,21 

0,12 

463.722,13 

28.709,03 

395,03 

1,39 

263.567,67 

27.757,06 

-556,94 

-1,97 

295.180,93 

28.959,47 

645,47 

2,28 

446.676,83 

28.391 ,63 

77,63 

2,74 

199.017,92 

27.789,84 

-524,16 

-1,85 

207.530,44 

28.130,67 

-183,33 

'-0,65 

121 .222,82 

27.810,09 

-503,91 

-1,78 

458.580,67 

28.141,33 

-172,67 

-0,61 

138.704,39 

28.313,02 

400,96 

1,40 


284.932,95 

- 

- 

- 

1,88 

- 

- 

- 

- 

0,01n.s. 

- 

- 

a) 

533,79 

1,90 

- 
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0 suntario dos resultados obtidos para os cinco pianos 

testados estao contidos na Tabela 4.15. 

* 

TABELA 4.15 

SUMflRIO DOS RESULTADOS ALCANgADOS PARA OS 5 PLANOS AMOSTRAIS TESTADOS 


S1STEKA DE AMOSTRACEM 

MEDIA 

(ha) 

DESVIO 

PADRKO 

(ha) 

C0EFICIENTE 
DE VAR!A?A0 
(*) 

PIFERENQA 

DE 

28.314 (ha) 

R.M.Q.E 

(ha) 

TAHAUHO 00 
SEG1ZNT0 
(km*} 

H9 DE UfllDADE 
POR AMOJTRA 
00 

PERCENTAGEM 
AMO 5 TRADA 
■(*) 

MEDIA 

{«) 

MAXIMA 

(ha) 

10 

10 

13,90 

28,500,10 

202,59 

2,92 

2,30 

-2065.51 

832,95 

20 

0 

22,22 

28.301,78 


3,47 


1882,82 

957,34 

30 

e 

33,33 

. 28.502,71 

491,17 

1,72 

1,40 

-1045,82 

526,20 

AO 

6 

33,33 

28.371 ,G8 

833,33 

SB 

2,50 


834,33 

CO 

c 

SO ,00 

28,313,02 


mm 

1,40 



1477,36 

533,79 


De urn modo geral , a utilizagao dos esquemas amostrais 
se'l ecionados forneceram estimativas exatas e precisas. Isto pode ser 
observado pelcs baixos valores das medias das diferongas relativas, 
ou da RaTz Media Quadrativa dos Erros (R.M.R.E.), e pelos baixos coe_ 
ficiente de variagao obtidos das 20 estimativas das areas tritTcolas 
para cada amostra selecionada. 

Da analise da selegao do tamanho do segmento e da de_ 
terminagao do numero de unidades por amostra, para os cinco pianos 
amostrais testados, o melbor desenipenho alcangado pel a utilizagao do 
metodo de estimativa pel a regressao foi obtido com o uso do segmento 
de 10 km?,:pois, alem de ter fornecido a estimativa da area ocupada 
com trigo tao precisa e exata quanto aos demais tamanhos de segmentos 
tais estimativas foram obtidas com a menor percentagem de area amos^ 
trada, 13,90%. 


Entre os menores tamanhos de amostras determinados pa_ 
ra os cinco pianos amostrais estudados,sel ecionou-se como melhor si£ 
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tema a amostra composta de 10 unidades com tamanho de segmentos de 
10 km 2 . 


Para verificar o desempenho do metodo de estimativa pje 
la regressao, procurou-se comparer os resultados obtidos por este me 
todo com aqueles obtidos pelo metodo de expansao direta, usando as me_s 
mas unidades amostrais. Na Tabela 4.16 estao conti dps os dados compara_ 
tivos dos metodos de estimativa pela regressao e de expansao direta. 

Conforme pode ser observado pel os resultados conti dos 
na Tablea 4.16, o metodo de estimativa pela regressao mostrou-se subs_ 
tancialmente mais eficiente quando comparado ao metodo de expansao di 
reta.Isto e demonstrado pela grandeza da eficiencia relative que va_ 
riou de 5,75 a 54,75. Em outras palavras, pode-se dizer que quando se 
utilizou dados auxiliares do MSS do LANDSAT para estimar a area da cul_ 
tura do trigo, atraves da estimativa pela regressao, o ganho em precj_ 
sao foi no mini mo de 5,75 ate 54,75 vezes maior do que aquele obtido 
pela expansao direta (utilizando somente fotografias aereas). A ca^ 
sa de variagao da Eficiencia Relativa (ER) esta fortemente ligada ao 
valor do coeficiene de determinagao (r 2 ) estimado para cada amostra, 
ou seja, quanto maior o r 2 estimado, menor a variancia estimada pelo 
metodo de regressao e, consequentemente, maior o valor da ER (Figura 
4.11). 


Por outro lado, as estimativas obtidas pela regressao 
tambem foram mais exatas, o que e demonstrado pela variagao das esti_ 
mativas das areas tritTcolas em ambos os metodos. Enquanto as estima 
tivas das areas triticolas obtidas pela regressao variaram de 
26.228,49 a 29.610 ha, estas mesmas estimativas obtidas pela expansao 
direta variaram de 19.370,74 a 35.513,21 ha. 
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lA 4.16 


'IVAS DAS AREAS DE TRIGO. QBTIDAS PELOS 


EXPANSE OIRETA 

EFICIENCIA 
RELATIVA (EP) 

T tt) (ha) 

Wed) 

23.946,52 

12.904.191,00 

13,68 

33.648,19 

4.919.018,99 

13,63 

24.407,93 

17.059.428,96 

54,75 

33.109,49 

9.807.997,74 

22,42 

27.779,83 

12.874.165,09 

12,69 

30.131,42 

12.123.997,82 

15,14 

30.497,47 ■ 

10.130.649,71 

17, eg 

35.513,21 

5.246,-, ,76 

8,71 

2G .240,83 

10.491 . 1 . 81 ,01 

24,19 

25,968,82 

13,198.144,85 

30,05 

32,910,91 

16.708.312,54 

16,53 

22.045,18 

7.850.574,42 

6,68 

33.382,01 

10.422.605,23 

5,82 

28.273,46 

8.969.806,35 

5,75 

24.835,82 

12.458.182,29 

50,59 

31.657,46 

11 .572. 221 , 63 

9,94 

19.370,74 

11 .302. 107, 97 

15,59 

27.626,47 

9.951 .465,25 

15,33 

27.490,19 

14.597.391,10 

14,25 • 

30.812,68 

23.269.900,02 

22,27 










EFICIENCIA RELATIVA { Eft} 
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Fig. 4.11 - Grafico que mostra o acrescimo da eficiencia relative com 
0 aumento do coeficiente de determinagao. 




CAPlTULO 5 


CONCLUSOES E RECOMENDAgOES 


De acordo com os resultados desta pesquisa, em que forani 
utilizados dados do MSS do LANDSAT como variavel auxiliar para estimar 
a area ocupada com trigo, chegou-se as seguintes conclusoes: 

1} 0 metodo de estimativa pela regressao foi selecionado em cons^ 
quencia de dois fatores: a) boa correlagao entre os dados obti. 
dos de aeronave e dados do MSS do LANDSAT; e b) a intersegao da 
reta que representou a relagao entre estes dados deu-se em pon. 
to diferente da origem (0,0); 

2) Para obter desempenho de estimativas semelhantes, a medida que 
se diminuiu o tamanho do segmento, a percentagem da area amos. 
trada tambem diminuiu; 

3) A menor percentagem de area amostrada que estimou a area de tri_ 
go com boa precisao e exatidao, atraves da estimativa pela re. 
gressao, foi de 13,90%, util izando-se como unidade amostral o 
segmento de 10 km 2 com 10 unidades por amostra. Com isto, houve 
uma economia de tempo e custo despondidos em cerca de 86% duraji 
te a fase de fctointerpretagao dos dados da aeronave; 

4) Com a utilizagao deste sistema amostral , pode-se obter estimati. 
vas por intervalo, o que e estatiscamente mais recomendavel do 
que a estimativa por ponto. Alem deste fator positive, este sis. 
tenia amostral forneceu estimativas com exatidao ^ 92,47% e com 
coeficiente de variagao das referidas estimativas menor que 5%. 

5) A area tritTcola pode ser obtida atraves de sistema de amostra^ 
gem sem utilizar dados do MSS do LANDSAT como variavel auxiliar, 
isto e, atraves da expansao direta (somente dados de aeronave) 
as estimativas foram menos exatas e precisas quando comparadas 
aquelas obtidas atraves do metodos de estimativa pela regressao. 
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6) No presente estudo, as areas obtidas pela classificagao digital, 
usando dados do MSS do LANDSAT, foram superestimadas para a maio 
ria dos segmentos e altamente correlacionadas com as mesmas areas 
obtidas atraves da interpretagao visual de fotografias aereas, iji 
dependentes do tamanho do segmento e do sentido considerado, po_ 
rem, em alguns segmentos ocorreu uma subestimativa. Os temas que 
mais se confudiram com o trigo foram mata e pastagem. 

Dada a experiencia adquirida durante este trabalho, reco. 

menda-se: 

1} para grandes areas onde as condigoes de relevo, tipos de solo, 

clima, epoca de plantio, etc. sao heterogeneas , que se faga uma 

estratificagao da area de interesse em subareas homogeneas(estra 
tos) e para cada estrado seja aplicada a estimativa pela regres^ 
sao, onde a area total sera obtida pela agregagao da area de ca_ 
da estrato. Desta forma, pode-se obter melhor precisao da estima 
tiva de area com o mesmo numero de umidades amostradas para esti_ 
mar a area sem a estratificagao; 

2) com relagao ao uso de fotografias aereas, que se teste uma alter 

nativa para substituigao de fotos aereas adquiridas no ano agrT 

cola da cultura de interesse, por coleta de informagoes direta_ 
mente no campo, com base em fotos aereas de anos anteriores, v_i_ 
sando com isto baixar o custo da aquisigao dos dados terrestres. 
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TABELA A.l 

» i I i i ■ Mii.« ■■ i.i.d 

DADOS UTILTZADOS NO CALCULO DAS AREAS TRITTCOLAS (FOTO AEREA E LANDSAT) 

PARA OS SEGHENTOS DE 10km 2 . 


N? DO 
SEGMENTO 

DENS I DADE DE 

CULTURA (DC) 

m 

AREA (ha) 

CLASSIFICADA 

CORRETAMENTE 

(1) 

DE NAO-TRIGO 
CLASS IF. COMO 
TRJGO (2) 

DE TRIGO NAO 
CLASS IF. COMO 
TRIGO (3) 

1 

44,25 

393,26 

158,24 

77,63 

2 

43,63 

397,15 

171,38 

67,18 

3 

46,78 

428,50 

150,78 

69,27 

4 

42,57 

358,31 

132,86 

94,66 

5 

60,61 

584,97 

1 88»6 7 

60,01 

6 

30,59 

292,59 

110,17 

32,84 

7 

53,85 

504,05 

132,86 

68,97 

8 

52,59 

511,82 

166,00 

47,77 

9 

39,42 

340,12 

180,04 

79,42 

10 

43,80 

417,16 

176,75 

48,96 

11 

55,50 

546,15 

100,92 

44,49 

12 

29,94 

249,63 

79,72 

68,97 

13 

, 64,48 

616,92 

74,04 

72,25 

14 

28,56 

267,55 

119,43 

38,22 

15 

63,67 

606,1 > 

105,07 

71,36 

16 

48,26 

438,36 

149,28 

75,24 

17 

56,06 

536,29 

123,31 

60,31 

18 

39,23 

379,22 

69,27 

38,22 

19 

17,17 

174,99 

95,84 

7,76 

20 

9,97 

93,17 

112,86 

10,15 

21 

26,38 

219,48 

81,21 

61,21 

22 

25,00 

243,66 

109,28 

22,39 

23 

29,24 

269,05 

72,25 

42,09 

24 

37,26 

363,41 

145,70 

31 ,17 


(continua) 
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Tabela A.l - Conti nuagao 



DENSIDADE DE 
CULTURA /OC) 
(*) 

AREA (ha) 

CLASS IFICADA 
CORRETAMENTE 
(1) 

DF. NffO-TRIGO 
CLASS IF. COHO 
TRIGO (2) 

DE TRIGO NAO 
CLASSIF. COMO 
TRIGO (3) 

25 

41,17 

406,99 

31,65 

42,39 

26 

19,72 

171,11 

25,67 

38,81 

27 

28,25 

258,29 

92,56 

42,40 

28 

45,68 

432,38 

83,30 

53,74 

29 

15,40 

151,99 

28,96 

11,94 

30 

36,08 

355,04 

54,64 

28,96 

31 

26,88 

245,07 

60,10 

40,99 

32 

63,59 

631,74 

40,14 

44,92 

33 

40,04 

349,37 

26,57 

76,73 

34 

18,60 

161,85 

32,84 

36,13 

35 

24,66 

186,22 

11,05 

77,03 

36 

0,64 

3,29 

18,21 

3,58 

37 

40,94 

380,15 

80,29 

55,50 

38 

43,32 

426,11 

36,72 

34,03 

39 

53,62 

546,29 

60,47 

24,33 

40 

45,72 

402,82 

123,01 

83,60 

41 

53,42 

538,38 

65,98 

30,16 

42 

37,29 

350,27 

44,78 

46,58 

43 

61,11 

585,74 

130,30 

64,61 

44 

60,52 

596,31 

37,92 

47,77 

45 

51,83 

514,79 

60,01 

36,72 

46 

45,79 

428,80 

49 86 

58,52 

47 

37,85 

369,68 

100,02 

33,14 

48 

31,54 

246,19 

40,17 

89,44 

49 

46,86 

472,40 

166,90 

26,26 

50 

32,02 

289,65 

78,82 

51,05 

51 

43,89 

401,33 

65,98 

65,68 

52 

18,80 

149,31 

59,41 

50,75 


( conti nua) 


c - x 
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Tabela A.l - Conclusao 


DO 

GMENTO 

DENSIDADE DE 
CULTURA (DC) 
(*) 

AREA (ha) 

CLASSIFICADA 

CORRETAMENTE 

0) 

DE NAO-TRIGO 
CLASS IF. COMO 
TRIGO (2) 

DE TRIGO NAO 
CLASS IF. COMO 
TRIGO (3) 

53 

52,02 

413,87 

30,45 

139,73 

54 

42,56 

320,62 

95,63 

132,36 

55 

50,96 

494,94 

70,52 

47,32 

56 

40,09 

391 ,18 

87,18 

35,53 

57 

22,17 

192,01 

70,31 

43,88 

58 

11,31 

84,81 

68,37 

35,53 

59 

22,48 

183,49 

30,31 

55,69 

60 

31,76 

227,54 

66,58 

110,47 

61 

53,62 

522,25 

65,39 

48,37 

62 

51,77 

464,04 

89,27 

86,88 

63 

41,47 

384,01 , 

59,41 

57,32 

64 

34,32 

320,11 

72,85 

45,08 

65 

6,17 

54,05 

39,41 

11,64 

66 

45,79 

406,70 

75,84 

80,62 

67 

25,23 

241,28 

40,90 

27,16 

68 

3,28 

30,46 

18,21 

4,48 

69 

26,18 

210,81 

49,86 

67,77 

70 

18,26 

157,67 

55,23 

36,72 

71 

1,93 

14,93 

29,26 

5,67 

72 

25,65 

246,95 

56,13 

26,02 
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TAB EL A A. 2 


DADOS UTILIZADOS NO CALCULO DAS AREAS TRITTCOLAS (FOTO AEREA E LANDSAT) 

PARA OS SEGMENTOS DE 20km 2 . 





AREA fhal 


N? DO 

SEGMENTO 

DENS I DADE DE 

CULTURA (DC) 
(*) 

: CLASSIFICADA 
CORRETAMENTE 
0) 

DE NAO-TRIGO 
CLASS IF. COMO 
TRIGO (2) 

DE TRIGO NAO 
CLASS IF. COMO 
TRIGO (3) 

1 

43,94 

790,43 

329,62 

144,81 

2 

44,67 

792,81 

283,64 

157,93 

3 

45,60 

877,61 

198,84 

92,85 

4 

41,98 

806,26 

177,01 

89,53 

5 

49,67 

949,11 

183,48 

107,93 

6 

45,36 

888,65 

110,73 

76,74 

7 

53,21 

1015,87 

298,86 

116,74 

8 

41 ,61 

757,29 

356,79 

128,38 

9 

42,72 

795,78 

180,64 

113,46 

10 

60,80 

1182,04 

168,22 

112,38 

11 

48,81 

943,59 

109,87 

95,24 

12 

34,69 

615,86 

140,19 

122,58 

13 

46,52 

884,48 

193,47 

105,68 

14 

55,97 

1044,54 

254,37 

146,60 

15 

47,64 

915,52 

192,58 

98,56 

16 

39,43 

762,06 

245,72 

77,31 

17 

31,24 

550,64 

125,39 

116,43 

18 

47,29 

734,48 

126,08 

272,09 

19 

13,69 

268,16 

208,70 

17,91 

20 

25,69 

463,14 

190,49 

83,60 

21 

33,25 

632,46 

217,95 

75,23 

22 

45,52 

885,90 

157,70 

82,85 

23 

16,73 

288,92 

138,68 

67,85 

24 

27,12 

411,02 

96,89 

166,16 
( conti nua) 
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Tabela A. 2 - Conclusao 

H AREA (ha) 

N? DENSIDADE DE CLASSIFICADA DE NAO-TRIGO DE TRIGO NAO 

SEGMENTO CULTURA (DC) CORRETAMENTE CLASSIF.COMO CLASSIF. COMO 



(*) 

(1) 

TRIGO (2) 

TRIGO (3) 

25 

30,28 

578,10 

57,32 

82,20 

26 

36,97 

690,67 

175,86 

96,14 

27 

25,75 

507,03 

83,60 

40,90 

28 

52,69 

986,29 

154,66 

135,25 

29 

37,89 

704,12 

132,26 

102,40 

30 

25,98 

460,75 

115,25 

92,26 

31 

45,23 

876,78 

100,24 

85,91 

32 

29,32 

511,21 

59,41 

112,86 

33 

12,65 

189,91 

29,26 

80,61 

34 

14,25 

271,74 

59,11 

31,64 

35 

22,22 

368,48 

105,09 

104,49 

36 

13,82 

261 ,88 

85,39 

31,69 
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TABELA A. 3 

DADOS UTILIZADOS NO CALCULO DAS AREAS TRIHCOLAS (FOTO AEREA E LANDSAT) 


PARA OS SEGMENTOS DE 30km 2 . 


BREA (ha) 


SEGMENTO 


DENSIDADE 

CULTURA 


44,89 

44.59 

45.94 
45,48 
48,62 
43,08 
57,82 

38.39 
52,24 

47.85 
40,92 

37.39 

17.75 

30.50 
37,74 

21 .85 

29.59 

32.39 

48.95 

28. 76 

43.50 
14,64 
18,23 

15,30 


DE CLASSIFICADA 
(DC) CORRETAMENTE 
1 ) 


1218,93 

1241,92 

1352,54 

1291,48 

1356.00 
1212,95 

1696.84 
1044,66 
1490,65 
1353,89 
1163,38 

883,80 

487,64 

876,13 

1090.01 

495.84 

836.40 

939.41 
1370,30 

780,86 

1226,15 

351,76 

482,56 

419,55 



480,40 

331 .70 

177.48 
233,77 

478.90 
357,39 

228.23 
190,05 
298,56 
341,86 

311.70 

185.49 

289.91 

327.23 
228,01 
165,26 
149,88 
166,90 
214,07 
188,10 
126,81 

62,10 

108,97 

140,62 


DE TRIGO NBO 
CLASSIF. COMO 
TRIGO (3 


214,08 

181.57 
113,86 

160.34 
196,16 

162,42 

149.10 

181.10 
177,03 
173,77 
142,99 

322.54 
79,12 
97,62 

114,84 

201,69 

108,37 

94,64 

192.57 

137.34 

162,64 

115.55 
99,41 

68,96 
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TABELA A. 4 


DADOS UTILIZADOS NO CALCULO DAS AREAS TRITlCOLAS fFOTO AEREA E LANDSAT 


>ARA OS SEGMENTOS DE 40km 2 . 


SEGMENTO 


DENSIDADE DE 
CULTURA (DC) 
(%) 

AREA (ha) 

CLASSIFICADA 

CORRETAMENTE 

(1) 

DE NAO-TRIGO 

CLASSIF.COMO 
TRIGO (2) 

DE TRIGO NAO 

CLASSIF. COMO 
TRIGO (3) 

48,58 

1806,30 

628,48 

261,55 

29,98 

1152,64 

402,17 

133,18 

37,84 

1454,88 

157,56 

155,89 

43,14 

1550,09 

640,43 

286,30 

40,83 

1507,67 

444,86 

230,18 

33,15 

1201,88 

235,27 

209,00 

44,15 

1673,38 

379,48 

206,31 

40,45 

1547,93 

410,53 

173,76 

19,20 

696,94 

112,86 

120,51 

51,45 

1988,31 

285,23 

201,91 

42,48 

1647,96 

403,42 

160,36 

33,48 

1258,03 

213,77 

166,89 

. 49,24 

1892,70 

293,35 

203,17 

24,04 

839,55 

264,07 

183,75 

30,06 

1072,60 

237,35 

206,90 

40,03 

1504,52 

250,95 

199,29 

37,20 

1145,52 

222,97 

438,25 

19,90 

722,63 

200,64 

124,50 
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TAB EL A A. 5 


DADOS UTILIZADOS NO CALCULO DAS AREAS TRITtCOLAS (FOTO AEREA E LANDSAT) 


ARA OS SEGMENTOS DE 60km 2 . 


46,75 

43,84 

51,88 

41,94 

34,99 

39,18 

39,33 

29,82 

36,55 

23,52 

33,59 

22,03 


Area fhai 

CLASS IFICADA 

DE .NAO-TRIGO 

DE TRIGO NAD 

corretamente 

CLASSIF.COMO 

CLASS IF. COMO 

CD 

TRIGO (2) 

TRIGO £3) 

2574,93 

959,30 

410,17 

2454,87 

689,09 

343,91 

3049,39 

405,71 

262,98 

2336,13 

423,82 

341,44 

1978,29 

588,47 

256,13 

2230,02 

699,09 

271,35 

2253,39 

539,71 

257,84 

1379,63 

350,75 

524,53 

2062,54 

276,69 

271,02 

1291,17 

229,00 

210,79 

1852,86 

323,04 

291 ,98 

1200,41 

328,72 

206,29 








APENDICE B 

PROGRAMA (EM LINGUAGEM FORTRAN) UTILIZADO PARA ESTIMAR A AREA DE TRIGO 

ATRAVES DO METODO DE REGRESSftO 
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PROGRAMA (EH LINGUAGEM FORTRAM) UTILIZADO PARA ESTIMAR A AREA DE TRIGQ 
ATRAVES DO METODO DE ESTIMATIVA PELA REGRESSflO 


FILE 

1 (KIND 

= REMOTE) 

FILE 

2(KIND 

= REMOTE) 

FILE 

3 ( KIND 

= PRINTER) 

FILE 

5 ( KIND 

= DISK,FILETYPE=7,TITLE="AUTO 

FILE 

6 (KIND 

= DISK,FILETYPE=7,TITLE="F0T0 


C - OBOETIVOS DO PROGRAMA: 

C - EXTRAIR AMOSTRAS ALEATORIAS SIMPLES SEM REPOS IQAO DE UMA POPULAQAO 

C - FITA. UTILIZA PARA EXTRACAO DAS AMOSTRAS A FUNQAO RANDOM. 

C - UTILIZA A ESTIMATIVA PELA REGRESSAO, CALCULA A AREA TOTAL E A VA 

C - RIANCIA DA ESTIMATIVA PARA UMA DETERMINADA CULTURA DENTRO DE UMA 

c - Area de interesse. 

C - XI - E UM NDMERO QUALQUER PARA INICIAR A FUNQAO RANDOM (i=l,n). 

C - VARIAVEIS: 

C - NEP NOMERO DE UNIDADES DA POPULAQAO 

C - IREP REPETigDES (1=1 ,n) 

C - NEA NOMERO DE ELEMENTOS DA AMOSTRA 

C - SA SOMA DAS AREAS CLASSIFICADAS PARA N SEGMENTOS AMOSTRAIS 

C - SF SOMA DAS AREAS INTERPRETADAS, VISUALMENTE DAS FOTOS PARA N 

C - SEGMENTOS AMOSTRADOS 

C - SDAF SOMATDRIA DOS DESVIOS DAS AREAS INTERPRETADAS VISUALMENTE VE 
C - ZES OS DESVIOS DAS AREAS CLASSIFICADAS DIGITALMENTE, OBTIDAS 

C - DAS UNIDADES AMOSTRADAS. 

C - SDA2 SOMA DOS DESVIOS AO QUADRADO PARA OS DADOS DIGITAIS DAS UNI_ 
C - DADES AMOSTRADAS. 

C - SDF2 SOMA DOS DESVIOS AO QUADRADO PARA OS DADOS DE FOTOS AEREAS 
C - DAS UNIDADES AMOSTRADAS. 

C - DMA2 DIFERENQA ENTRE A MEDIA POPULACIONAL, E AMOSTRAL, AO QUADRA 
DO DOS DADOS DE CLASSIFICAQflO DIGITAL. 
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C - AMEFOT 
C - 

C - AMGAUT 
C - 

C - YREG 
C - 

C - YTOT 
C - VAREST 


MEDIA DOS DADOS INTERPRETADOS VISUALMENTE PARA AS UN1DADES 
AMOSTRADAS 

MEDIA DOS DADOS CLASSIFICADOS DIGITALMENTE PARA AS UNIDA 
DES AMOSTRADAS 

AREA ME I DA ESTIMADA PELA REGRESSAO PARA AS UNIDADES AMOSTRA 
DAS 

AREA TOTAL ESTIMADA 
VARIANCIA ESTIMADA 


DIMENSION DE(300, ) ,DAM(100) , NAL ( 200) » IB(1 00 ) »DEF (300 ) 
WRITE(IJOOO) 

1000 F0RMAT(1 X, "QUAL 0 NOMERO DE ELEMENTOS DA POPULAQAO ?) 
READ ( 2 » /) NEP 
DO 2 1=1,300 

2 DE ( I ) =0 

DO 3 1=1,100 

3 DEF( I )=0 

READ(5,/) (DE(I) ,1=1 ,NEP) 

READ(6,/) (DEF(I) ,1=1 ,NEP) 

WRITE(1 ,1031 ) 

1031 F0RMAT( IX, "QUAL 0 NOMERO INICIAL.IZADOR? XI= ") 

READ(2,/)XI 

WRITE(1 ,1011 ) 

1011 FORMAT (IX, "QUANTAS REPETIQOES SERAO FEITAS? ") 
READ(2,/) IREP 
WRITE ( 3,1 0001 ) IREP 

10001 FORMAT( IX ,// ,1X , "NUM. DE REPETIQOES= ",13} 

DO 84 K=1 , IREP 
WRITE ( 1 ,1020) 

1020 FORMAT (IX, "QUANTOS ELEMENTOS TEM POR AMOSTRA?") 

READ ( 2 , / ) NEA 
DO 10 J=1 ,NEA 
A=RANDOM(XI) 

IA= A* NEP+0.5 
IF(IA.EQ.O) IA=1 
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I3(J)=IA 

XI=A*100Q 

CONTINUE. 

GO TO 21 

22 A=RADOM 

IA« A* NEp +0,5 
IF(IA.EQ.Q)IA*1 
IB(.J)=IA 

21 DO 23 0=1 , NEA-1 
DO 24 M=1 ,NEA-J 
IF( IB(0)EQ . IB(M+J) )G0 TO 22 
24 CONTINUE 

23 CONTINUE 
WRITE(3, 10003) 

10003 FORMAT(lx, "INDICES DAS AMOSTRAS",/) 
WRITE(3,/)CIBCJ)0=1,NEA) 

SA=0 

SF=0 

SDAF=0 

SDA2=0 

SDF2=0 

S=0 

DMA2=0 

DO 45 J=1 ,NEA 
SA=SA+DECIBC0).1 
SF=SF+DEF( IB toil 
44 CONTINUE 
AMEF0T=SF/FL0ATCNEA1 
AMEAUT=SA/FL0AT(NEA1 

C- CALCULO DA MEDIA GLOBAL DA CLASSIFICAgAO AUTOMATICA 

DO 55 J=1,NEP 

S=S+DECJ1 

55 CONTINUE 

MAGAUT=S/FLOATtNEPl 

DO 65 J=1 ,NEA 
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$DAF«SDAF+( (DE(IB(J J ) -AMEAUT)*(DEF(,IB(J ) ) -AMEF0T1) 

SDA2=SDA2+(_{DEf IB(Jl) -AMEAUT).*(OE(IB( J}) -AMEAUTJ.) 
SDF2 bSDF2+C(PEFCIBCJ)1-AMEF0T)*CDEF(IBCJ).)-AMEF0T}) 

65 CONTINUE 
B=SDAF/SDA2 

YREG=AMEFOT+0*CAMGAUT-AMEAUT) 

YTOT=NEP * YREG 

IFCNEA.EQ.21G0 TO 76 

DMA2= ( AMGAUT -AMEAUT ) *(AMGAUT- AM EAUT 1 

R 2 = (SDAF * SDAF) / ( SDA2 * SDF2) 

VAREST»(NEP*NEP7FL0AT{NEA)*(NEP-NEA)/NEP*SDF2*O~R 2 )/(NEA-2))* 

*(1+1./NEA+DMA2/SDA2) 

76 WRITE (3,1070)NEP,NEA,YT0T,VAREST,R 2 ,SDF2,SDA2,DMA2 
1070 F0RMAT(IX j 4(/),1X,"PLAN0M3,/,1X,"ELEMENT0SM3,// 1 1X, ,, AREA 
*ESTIMADA= ,i ,F 15.3,/,1X,"VARIANCIA ESTIMADA= ",F15.37,1X,"R 2 =V8.6,* 
*/>lX,"SDF2=" ,F1 5.3,/>TX,"SDA2 c",F 15 .3,/»l X, " DMA2= , Fl 5 . 3 ) 

84 CONTINUE 

STOP 

END 



